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Abstract 

Wc present a simple empirical parameterization of the x- and t-dcpcndcncc of generalized parton 
distributions at zero skewness, using forward parton distributions as input. A fit to experimental 
data for the Dirac, Pauli and axial form factors of the nucleon allows us to discuss quantitatively the 
interplay between longitudinal and transverse partonic degrees of freedom in the nucleon ( "nucleon 
tomography" ) . In particular we obtain the transverse distribution of valence quarks at given mo- 
mentum fraction x. We calculate various moments of the distributions, including the form factors 
that appear in the handbag approximation to wide-angle Compton scattering. This allows us to 
estimate the minimal momentum transfer required for reliable predictions in that approach to be 
around \t\ ~ 3GcV 2 . We also evaluate the valence contributions to the energy-momentum form 
factors entering Ji's sum rule. 
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1 Introduction 



In recent years hard exclusive reactions have found increased attention because of new theoretical 
developments. For a number of such reactions, for instance deeply virtual El or wide-angle [U E] 
Compton scattering off the nucleon, the scattering amplitudes factorize into partonic subprocesses 
and soft hadronic matrix elements, called generalized parton distributions (GPDs) HH Ej . While 
the partonic subprocesses can be evaluated in perturbation theory, calculation of GPDs requires 
non-perturbative methods. As an approach starting from first principles, lattice QCD is suitable 
for evaluating x-moments of GPDs. Interesting results on the first few moments have recently been 
obtained [8ll9"l ll0) lll). although such calculations still come with important systematic uncertainties. 

At present one is therefore unable to get along without models for the GPDs in order to describe 
data on hard exclusive reactions. These models are frequently simple parameterizations, constrained 
by general symmetry properties, by reduction formulas which express that certain GPDs become the 
usual parton densities in the forward limit of vanishing momentum transfer, and by the sum rules 
which express that the integrals of quark GPDs over x give the contributions of these quarks to 
the elastic form factors of the nucleon. Other attempts to model GPDs are for instance based on 
effective quark- model Lagrangians |121 113) or on the concept of constituent quarks |14) 115) . Another 
approach constructs GPDs from the partonic degrees of freedom through the overlap of light-cone 
wave functions [5} 116) ITT). For large x or for large momentum transfer, where essentially only the 
valence Fock state matters, this concept provides reasonable results as comparison with the measured 
parton densities and form factors reveals 0. Recent reviews of the general properties of GPDs and 
efforts to model them can be found in )18l 119) . 

A long-term goal is the (almost) model-independent extraction of GPDs from experimental data, 
in analogy to the determination of the usual parton densities from hard inclusive processes performed 
for instance in |2U1 121) . For GPDs this is admittedly a very challenging task since they are functions 
of three kinematic variables: the average momentum fraction x of the partons, the skewness £, and 
the invariant momentum transfer t. Furthermore, compared to inclusive processes, we will for quite 
some time have much less data on exclusive reactions at our disposal. In the near future we can 
aim at suitable parameterizations of the GPDs, with a few free parameters adjusted to data. An 
attempt of such a parameterization will be presented in this work. The proposed parameterization 
is physically motivated on one hand by Regge phenomenology in the limit x — > and, on the other 
hand, by the physical intuition gained in the impact parameter representation of GPDs. The free 
parameters of this ansatz are fitted to the experimental data on the Dirac and Pauli form factors 
of the nucleon, exploiting the sum rules mentioned above. Only the valence quark combinations of 
parton distributions are constrained by these observables, but the combination of proton and neutron 
form factors allows for a separation of u and d quark distributions in a certain kinematical range. 
Similarly, the GPDs for polarized quarks can be constrained by the axial-vector form factor. The 
kinematic dependence of GPDs involves two aspects: the interplay between the two independent 
longitudinal momentum fractions x and £, and the interplay between the longitudinal variables and 
t. As a consequence of Lorentz invariance, the ^-dependence of GPDs drops out in the sum rules for 
the form factors of the quark vector and axial- vector currents ^ . We therefore restrict our study to 
GPDs at £ = and concentrate on the correlation between the variables x and t. This can be done 
in a physically very intuitive representation: after a Fourier transform to impact parameter space, 
GPDs at zero skewness £ describe the joint density to find a parton at a given longitudinal momentum 
fraction x and a given transverse separation from the center of the nucleon in the infinite-momentum 
frame [22) . 

For the valence part of the GPDs H and H we will use their forward limits q(x) and Aq(x) as 
an input to our parameterization, and our fit to form factor data will allow us to extract the average 
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impact parameter (b 2 ) x of valence quarks with a given momentum fraction. For the proton- helicity 
dependent GPD E our task is more difficult since we need to make an ansatz for both its forward limit 
e(x) and for its dependence on t or on impact parameter. Lack of sufficient data on the pseudoscalar 
nucleon form factor prevents a similar study for the GPD E at present. 

With the GPDs at hand we are in the position to compute various moments and compare them for 
instance with lattice QCD results. We can in particular evaluate the contribution of valence quarks 
to Ji's angular momentum sum rule pQ. In impact parameter space our results can be turned into 
"tomographic" images of the proton as suggested in [221 |2H1 HU • We shall also discuss wide-angle 
Compton scattering in some detail. The soft hadronic matrix elements appearing in the soft handbag 
description of this process are new form factors, which are expressed as 1 /x-moments of GPDs jlj Ej 
and can also be evaluated from our results. Comparison of the corresponding observables with preci- 
sion data that have been taken at Jefferson Lab and will be published shortly will subsequently 
allow for an examination of our theoretical understanding of wide-angle Compton scattering. 

This paper is organized as follows: In Sect. El we recall some properties of GPD at zero skewness. 
In Sect. we present the physical motivation of our parameterization and analyze the GPD H. The 
corresponding analyses of the GPDs H and E are described in Sects. El and El respectively. Various 
properties of our GPDs are shown in Sect. El In Sect. [7| we discuss wide-angle Compton scattering in 
the soft handbag approach and investigate the corresponding form factors. The paper ends with our 
conclusions in Sect. |H1 In two appendices we provide details on the nucleon form factor data we have 
used ( App. lA|) and collect all fit results (App. Ib|) . 

2 Generalized parton distributions at zero skewness 

Let us start by recalling some properties of generalized parton distributions at £ = 0. We use Ji's 
definitions of GPDs and their arguments ^ and for simplicity suppress the argument £, writing 
H(x,t) instead of H(x,£ = 0,t), H(x,t) instead of H(x,£ = 0,t), etc. 
Let us first concentrate on the combination 

H*(x, t, /i 2 ) = H*(x, t, /i 2 ) + m{-x, t, /x 2 ) (1) 

of the quark helicity averaged GPDs for flavor q in the proton. This is the combination entering the 
proton and neutron Dirac form factors as 

Ffit) = f dx (§ftf (x, t, /i 2 ) - ii?„ d (x, t, fi 2 )) , 

F?(t) = ^dx^H^xX^-lKix,^^)). (2) 

J 

In the forward limit t = the distribution H%(x,0) becomes the valence quark density q v (x) = 
q(x) — q(x). In both Ff and F™ we have neglected the contribution e a H^ from strange quarks: the 
difference of strange and antistrange distributions in the nucleon is not large |26j and the strange 
contribution to nucleon form factors at low t is seen to be small in neutral-current elastic scattering 
|27j . In (J2) and @ we have displayed the dependence of the GPDs on the factorization scale /i, which 
we will often omit for ease of notation. We will also use the notation Ff(t) = Jq dx H$(x, t) for the 
individual quark flavor contributions to the Dirac form factor. 

As shown by Burkardt |22| I24j , a density interpretation of GPDs at £ = is obtained in the mixed 
representation of longitudinal momentum and transverse position in the infinite-momentum frame. 
In particular, 

qv(x,b) = J ^0L e - ibA H«(x,t = -A 2 ) (3) 
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gives the probability to find a quark with longitudinal momentum fraction x and impact parameter 
b minus the corresponding probability to find an antiquark, where we reserve boldface notation for 
two-dimensional vectors in the transverse plane. The average impact parameter of this distribution 
at given x is 

jd 2 bb 2 q v (x,b) _B 



(4) 

t=o 

Since q v (x,b) is a difference of probabilities, (b 2 )% is not an average in the strict sense. It gives 
however the typical value of b 2 in q v (x,b) as long as this distribution is positive (which is the case 
for the parameterizations we will use, and which is generally expected when x is sufficiently large to 
neglect antiquarks compared with quarks). 

GPDs can be written as the overlap of light-cone wave functions. In impact parameter space this 
representation has an especially simple form: 

q v (x,b) = ^47r)^- 1 /ndx i ff(l-2x i )/nd 2 6*^(E^ 

N,/3 i i i i 

i ~ 1 2 

^T/j5(x-X J -)<y (2) ( & - b j)|^JV/?( £C t: & t)| • ( 5 ) 
3 

The index i runs over the N partons in a given Fock state, whose quantum numbers are collectively 
denoted by the index /3, and * s the light-cone wave function of this Fock state in a proton 
with positive helicity. This impact-parameter wave function is obtained from the wave function in 
momentum space by a Fourier transform as given in |28j . The index j singles out the struck parton 
and runs over all quarks or antiquarks with flavor q, with rjj = 1 for quarks and rjj = — 1 for antiquarks. 

As explained in |24| the impact parameter b in q v (x, b) is the transverse distance between the struck 
parton and the center of momentum of the hadron (see Fig. y). The latter is the average transverse 
position of the partons in the hadron with weights given by the parton momentum fractions. It was 
chosen to be the origin in ©, so that the transverse positions 6j and momentum fractions Xi of the 
partons satisfy J2i x i^i = 0. The center of momentum of the spectator partons is easily identified as 
— bx/(l — x). The relative distance b/(l — x) between the struck parton and the spectator system 
provides an estimate of the size of the hadron as a whole, and we denote its average square by 



It does however not account for the spatial extension of the spectator system itself, which remains 
unaccessible in quantities like GPDs at zero skewness, where only one single parton within a hadron 
is probed. From Fig. ^ one readily sees that d q provides a lower limit on the transverse size of the 
hadron. This quantity has also been considered in recent work on color transparency |29l !3(Jj . 

Just as the usual quark densities, GPDs depend on the factorization scale at which the partons 
are resolved. For £ = the evolution in \i is described by the usual DGLAP equation, which for the 
valence combination H$ reads 

» 2 ^H«(x,t,v 2 ) = f x y [pQ] + fl?(M,/* Z ) , (7) 

where [ ]+ denotes the usual plus-distribution, and where the kernel reads P(z) = ^ Cp at 
leading order in the strong coupling. We note that the situation for fi 2 <^ —t is rather subtle; it 
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Figure 1: A three-quark configuration in the mixed representation of definite transverse position and 
definite plus-momentum, x denotes the plus-momentum fraction of the upper quark with respect to 
the nucleon. The dashed line indicates the center of momentum of the two lower quarks and the thick 
solid line the center of momentum of the proton. 

will be discussed in some detail in Sect. 17.21 Dividing (J7J) by H%(x,t) and subsequently taking the 
derivative in t at t = we obtain an evolution equation for the average impact parameter: 




where the plus-prescription is no longer needed since the term in square brackets vanishes at z = x. We 
see that the average impact parameter decreases with // for all x, provided that (fo 2 )| is a decreasing 
function of x. 

Notice that the right-hand sides of (J2J) must be evaluated at a particular resolution scale (M, whereas 
the left-hand sides are the form factors of a conserved current and hence independent of the scale 
fi where the current is renormalized. Physically speaking, the transverse distribution of quarks of 
a given momentum fraction x is modified by the parton splitting processes that underly DGLAP 
evolution. It hence depends on the spatial resolution at which quarks are probed, whereas the 
transverse distribution of charge described by the electromagnetic form factors does not |2H] • 

For the quark helicity dependent GPDs we define the valence combination 

H*(x, t, fi 2 ) = H*(x, t, fx 2 ) - H*(-x, t, fi 2 ) , (9) 

whose forward limit is Aq v {x) = Aq{x) — Aq(x). The impact parameter distribution Aq v (x, b) is then 
defined in analogy to (JHJ and again can be interpreted as a difference of probability densities in x and 
b space. It has a wave function representation akin to ©, with an extra minus sign in front of the 
squared wave function if the struck quark or antiquark has negative helicity. The evolution of H$ in 
the scale fi is described by a DGLAP equation as in (|7|). with an evolution kernel that is identical 
to the one for H$ at leading order in a s . The properties of the proton helicity flip GPD E will be 
discussed in Sect. |5J 

3 A parameterization for the unpolarized GPD H v (x,t) 
3.1 Physical motivation 

We now develop a parameterization for H^(x,t) and H^(x,t). For its functional form we will use 
theoretical guidance in the regions of very small and very large x and then attempt a suitable inter- 
polation for intermediate x. We will fit this parameterization to the data on the nucleon Dirac form 
factors Ff(t) and Ff(i). 
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For small t and very small x one can expect Regge behavior of H v (x,t), employing the same 
argument as in the well-known case of forward parton densities |31| . The simplest form of Regge 
behavior is the dominance of a single Regge pole, 



H v (x,t) ~ c(t) (-fy 



XQ\a(t) 
X 



X \a(0) 



— exp 



a' log — + bo) t 

x / 



(10) 



where in the second step we have used a simple parameterization c(t) = exp(6o t) for the small-i region 
and a linear form of the Regge trajectory a(t) = a(0) + a't. The leading Regge trajectories with the 
quantum numbers of H v are those of the oj and the p. They can be phenomenologically determined 
from suitable hadronic cross sections and from the Chew-Frautschi plot (showing the spin of a meson 
versus its squared mass, which we take from |32| ) . From the masses of u;(782) and ^3(1670) one obtains 
a linear trajectory a u (t) = 0.42 + 0.95 GeV -2 i, whose intercept at t = agrees well with the value 
extracted from o- tot (pp) — crtotipp) in I3H] The masses of /?(770) and p3(1690) give a linear trajectory 
a p (t) = 0.48 + 0.88 GeV -2 1, in good agreement with the intercept from CFtot{^~p) — o"tot(^ + p) and 
with the trajectory extracted from the data on da/dt(ir~p — > 7r°n) up to about \t\ ~ 0.3 GeV 2 |34j . 

We emphasize that (|10[) is not a prediction of Regge theory, but rather corresponds to a simple 
form of Regge phenomenology: on one hand one expects subleading Regge trajectories to become 
important if x is not sufficiently small, and on the other hand the importance of Regge cuts, which 
lead to a more complicated behavior on x and t, is notoriously difficult to determine without further 
assumptions. To assess how well the ansatz IJlOjl fares at t = we have investigated the CTEQ6M 
parton densities [2U] at fj, = 2 GeV and found that for 10~ 5 < x < 10 -2 one has u v + d v ? 



-0.435 



(u v -d v )/{u v + d v 



-0.07 



and 

both within 1% accuracy. Scanning the 40 sets of parton densities given 
by CTEQ as error estimates, we found an exponent in the power-law for u v + d v between —0.38 and 
—0.495 and a corresponding exponent for (u v — d v )/(u v + d v ) between —0.165 and 0.13. Similar values 
are found when taking the distributions at scales fi = 1 GeV or fi = 4 GeV. We conclude that a simple 
Regge pole ansatz with an intercept taken from the phenomenology of soft hadronic interactions is in 
fair agreement with valence quark distributions at low factorization scale, and assume in the following 
that this description generalizes to small finite —t. Note that the form fjlOjl translates into an average 
impact parameter (b 2 ) x diverging like log(l/x) at small x. A physical mechanism that gives such 
a behavior is Gribov diffusion, the generation of small- :c partons through a cascade of branching 
processes |33j . 

As x increases, the struck parton takes more and more weight in the center of momentum J2i x %^% 
of all partons, so that the distribution in b should become more and more narrow (HHl- This means 
that the t-dependence of GPDs should become less steep with increasing x. If the average distance 
d q {x) between the struck quark and the center of momentum of the spectators is to remain finite, 
which one may expect for a system subject to confinement, then (b 2 ) x must vanish at least like (1 — x) 2 
in the limit x — > 1 |HZI- The actual limiting behavior of (b 2 ) x in QCD remains unknown. Certainly the 
impact parameter dependence of GPDs at large x contains interesting information about the dynamics 
of confinement, and we shall see how much information on this dependence can be extracted from 
existing form factor data. 

In our parameterization we will make an exponential ansatz for the t-dependence: 



H${x, t) = q v (x) exp [tfJx)] . 
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The function f q (x) parameterizes how the profile of the quark distribution in the impact parameter 
plane changes with X ? ELS IS readily seen from 



q v (x,b) 



4^ fq(x) 



4/«(a0 



{b 2 )% = 4f q (x) . 



(12) 
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Apart from being suitable for analytic calculations, an exponential t-dependence of H%(x, t) guarantees 
that q v (x, b) is positive. It ensures a rapid falloff at any x as —t becomes large, and it readily matches 
with the Regge form QlUj) for small x and —t if we impose 



f q (x) -> a' log - + B q 

x 



for x — > 



(13) 



The t-slope at x = xq is obtained as 60 = a ' l°g(l/a;o) + B q . I n our fits we will explore a possible 
flavor dependence of B q but keep a' flavor independent, as suggested by Regge phenomenology. We 
are aware that our exponential ansatz (jllj) has no rigorous theoretical backing, and we shall explore 
alternative forms of the t-dependence in Sect. 13.61 Let us emphasize already here that we cannot 
trust the detailed form of our extracted GPDs in the region —t S> l/f q (x), where according to 
they are exponentially small and thus give only a tiny contribution to the form factor integrals (J2J. 
In particular we claim no validity of our ansatz at —t of several GeV 2 and very small x, where its 
motivation from Regge phenomenology does indeed not apply. 

For large x one can expect that the overlap representation (JSJ is dominated by Fock states with 
few partons. In [S] we have evaluated GPDs from model wave functions for the lowest Fock states, 
whose dependence on the transverse parton momenta or impact parameters was taken as Gaussian, 



tp oc exp 



k 2 



ip oc exp 



1 



4^2 E -A 2 



(14) 



a form going back to [3H1 and explored in detail for the nucleon in |39j . With a parameter a 2 = 
0.72 GeV -2 or somewhat larger, this model allowed a fair description of unpolarized and polarized u 
and d quark densities for x>0.6 to 0.7 and of Ff(t) for — i>4GeV 2 . The resulting GPDs at £ = 



take the form given in (|lljl and (|12|) with (6 



2a 2 (l - x)/x. The average distance between the 



struck quark and the spectators hence diverges like d q (x) ~ (1 — x)~ x l 2 in the limit x — » 1. Indeed 
the impact parameter form of the model wave functions (|14|) allows b 2 to grow like l/x{ when the 
spectators become soft. In the limit x — ► 1 such a behavior is difficult to reconcile with confinement 
as was pointed out in [37] , and one aim of our study here is to explore quantitatively at which x the 
behavior of such wave functions becomes physically suspect. In our ansatz 1)11)1 for the valence GPDs 
we will impose the constraint 



f q {x) - A q (l 



for x — ► 1 



(15) 



either with n = 1 as in the model just discussed, or with n = 2, which corresponds to d q (x) tending 
to a constant at x = 1. 

The ansatz (|11|) must be made at a particular factorization scale /i and may work better for some 
scales than for others. Let us see that the limiting behavior we take for (b 2 ) x at small and at large x 
retains its form under leading-order DGLAP evolution. To be more precise, let us first assume that 
(b 2 ) x 4a'log(l/x) + 4B and q v (x) « cx~ a with a > at small x and for a given /i. We need not 
take the mathematical limit of x — > but only require these forms to be good approximations in a 
range of x where the small-x approximations of the following arguments are numerically adequate. 
With the evolution equation (jHJ) for (b 2 ) x and the leading-order evolution kernel we have 



2ir 





fX\ 2 


[ dz 


Wl) 


IX 





Qv(z) {b 2 
q v (x) 



(16) 



Let 5 be a fixed value of x below which q v (x) and (b 2 ) x can be approximated as stated above. For 
z > 5 the integrand behaves like x a logx for x — > 0, so that the integral over z from 5 to 1 gives a 
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vanishing contribution to the right-hand side. We can hence approximate 



o s Cf 



2tt j x 



dz 



(~ 

z 



\(J \og[zjx) 

Z' 



-4a 



, ol s Cf 
2vr 



/S/x 
du(l + u~ 



u 



z — X 
logu 



(17) 



which tends to a negative constant for x — > 0. The divergent part 4c/ log(l/x) of (b 2 ) x is hence [i 
independent, whereas the constant AB will decrease with fi. Our argument can be generalized to 
other forms of q v (x) at small x, for instance to a sum J2 k CkX~ ak with a^ > 0. 

Concerning the large- x behavior, one can show that a form (b 2 ) x ~ 4A(1 — x) n with n > is 
stable under leading-order DGLAP evolution, provided that the forward densities at a given [i behave 
as q v (x) ~ (1 — x)P . More precisely, the coefficient A decreases with /i, whereas the power n remains 
stable. To see this one starts with the evolution equation Qlfijl. replaces q v and (b 2 ) on the right-hand 
side with the approximations just given, and Taylor expands the evolution kernel to leading order in 
(1 — x). The result is 



-4A(1 



v 



7T 



1 I 

dv — 

l—V 



(18) 



The leading x-dependence is hence in (1 — x) n on both sides, and one obtains an equation for the 
evolution of its coefficient 4^4 with \i. Our finding is in line with a numerical study of pion GPDs 
by Vogt |4(Jj . who found that in a finite interval of large x the form (|llj) with j(x) = ^a 2 (l — x)jx 
is approximately stable under DGLAP evolution, with a moderate decrease with fx of the parameter 
a 2 (/x). For the Gaussian model wave functions giving this form of GPDs, a decrease of a 2 entails a 
decreasing probability of the corresponding lowest Fock states. This is in agreement with physical 
intuition: at higher resolution scale \x one resolves more and more partons in the hadron, and to find 
a configuration with only a few partons becomes less likely. 

The exponential t-dependence tf 1 If) of GPDs is generally not stable under DGLAP evolution. To 
see this let us consider the Taylor expansion 



logH%(x,t) =logq v (x) +t 



0_ 

Of 



log H*(x,t) 



J t=o 



1 2 



dt 2 



logfl?(M) 



+ ... , 



(19) 



J t=o 



which ends after the linear term if the t-dependence of H% is exponential. Dividing (J7J) by H%(x,t) 
and taking the second derivative in t we obtain the scale dependence of the quadratic term in (|19l) as 



2 1 
d/i 2 



^ ,.,2 



dt 2 



log^(M) 



1 



J i=0 



1 ^ z p /: 



q v (x) Jx z \z 



lv(z) 



(-log H%(z,t) - -log H$(x,tj) + w (log H%(z,t) - log H^(x,t] 



(20) 



t=o 



If at a given scale /x the t-dependence is exponential, then the right-hand side of (|2U|) is positive so 
that the quadratic term in ()19|) becomes positive as one evolves to a higher scale. 

In the small- x limit we do however find approximate stability under evolution. With the same 
assumptions and approximations that led to (|17|) we get 



2 d 



9!. 

at 2 



logfl^(x) 



(a ) 



2tt 



/d/a; 
du(l + n~ 2 )ii 



U - 1 



(21) 
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if the second derivative in t of \ogH 2 vanishes at the scale \i. The change with fj, of the quadratic 
term in the Taylor expansion (|19|) is then of order (ta 1 ) 2 . For moderate values of tot' this is small 
compared with the linear term to'log(l/x). 

In the large- x limit we get in analogy to 1)18(1 

^ 2 (l-x) 2 "^ fdvvP {l ~ Vn? (22) 

if at the scale \i the t-dependence of H q is exponential. The change with \i of the quadratic term in 
the Taylor expansion (J19|) is then parametrically of order {tA) 2 (\ — x) 2n . This is not small compared 
with the linear term tA(l — x) n if the latter is of order 1. Numerically however the u-integral in (|22|) 
is rather small, namely 0.05 for n = 1 and 0.14 for n = 2 if /3 = 3. For large x we thus find that the 
departure from the exponential behavior (|11|) under evolution should not be too strong in the t-region 
where the exponent tf(x) does not take too large negative values. This is again in agreement with 
the numerical study of Vogt [3U] mentioned above. 

So far we have considered the valence combination H v of quark GPDs. Let us briefly comment 
on what one would expect for the "sea quark" distributions H (x, t) at x < 0. At large x, the wave 
function overlap picture suggests a different impact parameter or t-dependence than for the valence 
distribution, because sea quark distributions require Fock states with at least one qq pair in addition 
to the minimal three-quark configuration. In the small- x limit one may expect a form as in (|10l) at 
low scale fi, given that the leading p, u, a\ and f\ Regge trajectories all have approximately the 
same a(t). The situation for sea quarks is however more complicated because the singlet combination 
^Z q [H q {x, t) — H q (—x, t)] mixes under DGLAP evolution with the gluon GPD H 9 (x, t), whose small-x 
behavior is dominated by Pomeron exchange. It is well-known that for the forward quark densities this 
leads to a drastic modification of the small-x behavior as one increases the scale \x even to moderate 
values of a few GeV, and one cannot exclude similarly strong modifications for the parameter a' of 
sea quarks. There is no data for form factors which might constrain the sea quark GPDs at £ = in a 
similar fashion as the electromagnetic form factors constrain H q . To investigate the sea quark sector 
one will rather rely on measurements of exclusive processes like deeply virtual Compton scattering or 
meson electroproduction at small t, where GPDs at nonzero £ are accessible. 

3.2 Selecting a profile function 

The assumption of an exponential t-dependence and the parameterization for the profile function 
f q (x) in 1)11(1 represent a source of theoretical bias, which translates into a systematic error in the 
determination of GPDs from experimental data. To gain a feeling for this error we will carefully 
compare different parameterizations. Our criteria for a good parameterization are: 

• simplicity, 

• consistency with theoretical and phenomenological constraints, 

• easy physical interpretation of parameters, if possible, 

• stability with respect to variations of the forward parton densities. 

In this section we discuss a few examples, including the default parameterization we will use in 
the remainder of the paper. For each parameterization we fix the free parameters by a x 2 fit to 
the experimental data on Ff and -F{\ For the forward densities q v {x) we will use the CTEQ6M 
distributions J^Uj at \x = 2 GeV as a default, where the choice of scale is a compromise between being 
large enough for q v (x) to be rather directly fixed by data and small enough to make contact with soft 
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Figure 2: Region of x (white region) which accounts for 90% of Ff(t) in the best fit to (|23j) at 
fi = 2GeV. The upper and lower shaded x-regions each account for 5% of Ff(t), see (|24j) . The thick 
line shows the average (x)t defined in (|25jl. 



physics like conventional Regge phenomenology. Tables with the results of our fits are collected in 
App. [Bl and details of our data selection and error treatment are given in App. 

The simplest form of the profile function f q {x) satisfying our constraints (|13[) and l)15J) with n = 1 
is actually f q (x) = a' log(l/x) itself. Such a Regge behavior of the GPDs has already been mentioned 
in an d [Ml an d was explored in some detail in ^Sj. One can however not expect this simple 
form, where one and the same parameter describes physics at small and large x, to work beyond a 
rough accuracy. Note that even in the small- x limit this form is special since it fixes the parameter 
bo in to be c/log(l/xo). As a simple extension of this ansatz one may try 

f q (x) = a'log- + (A q -a')(l-x), (23) 
x 

which is still in agreement with (|13|) and ()15|) . A very good fit (x 2 /d.o.f. = 1.34) of the nucleon Dirac 
form factors can then be obtained with three free parameters, a' , A u and Ad, see TableEJ The fitted 
value a' = (1.38±0.01) GeV~ 2 is however significantly larger than what Regge phenomenology would 
lead one to expect. This disagreement becomes even stronger if we take the forward parton densities 
at n = 1 GeV instead of 2GeV. The fit then gives a' = (1.65 ± 0.02) GeV" 2 with xVd.o.f. = 1.14. 

To better understand the situation we first determine the region of x in @ to which our fit is 
actually sensitive. To this end we consider the minimal and maximal value of x which is needed to 
account for 95% of the form factor in the sum rule, 

i 

J dxY,e q H q v (x,t) 

3-min(^) 
^max(^) 

f dxJ2e q H?,{x,t) 
o i 

where e u = | and = — |. We concentrate on the proton form factor for this purpose, which is the 
most important input to our fit given the available data. A related quantity is the average value of x 



95% • Ff(t) , 



(24) 
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Figure 3: The function l q (x) from (|26|) for the best fit to (|23|) with fj, = 2GeV. The contributions 
from terms going with a' (dashed line) and A q (dotted line) are shown separately. 



in the form factor integral, given by 

, . Io dx J2 q e q xH^(x,t) 

\ X /t = J g • \^) 

Jo dx Y, q e q H%(x,t) 

In Fig. [21 we show x m i n , x max and (x)t obtained for the best fit to (|23j) as a function of \t\. The 
relevant x-range moves towards higher values for increasing momentum transfer. With the existing 
data on Ff going up to -t = 31.2 GeV 2 , the region of x where the profile function can be constrained 
by our fit goes up to about 0.8 or 0.9. On the other extreme we have x m - m « 1.4 x 10 -3 at t = 0. 
Clearly, the form factor integrals are only sensitive to the small- x behavior of H% if t is also small, as 
we anticipated below 1)13(1 . 

Let us now see which of the parameters in our fit is most important in the profile function at 
given x. In Fig. we show the quantities 

1,00 = f M (26) 

resulting from our fit, where we have divided out the factor log(l/x) in order to have a finite quantity 
in the limit x — > 0. We also show the individual contributions to l q from the terms a' [log(l/x) — (1— x)\ 
and A q (l — x) in f q and see that the value of a' controls the profile function in almost the entire x range 
for u quarks and in a substantial region for d quarks. The fitted value of a' thus reflects dynamics at 
both small and large x (in the fit it has to find a compromise between these regions). We can hardly 
expect it to give a good representation of the physics in the region where Regge phenomenology is 
relevant, say for x < 10 -2 . 

The simplest profile function satisfying the constraints (|13|) and (|15|) with n = 2 is f q (x) = 
a'(l — x) log(l/x), which has been proposed in and used for numerical studies in |24| , An obvious 
extension of this ansatz is 

f q (x) = q'(1 - x) log l + (A q - c/)(l - xf . (27) 

A fit with free parameters a', A u and does not give a good description of the form factor data, 
see Table El Having an overall x 2 /d.o.f. = 5.25 it systematically overshoots the Ff data at \t\ above 
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Figure 4: The function l q {x) from ()26|) for the best fit to ()29|) with /i = 2GeV. The contributions 
from terms going with a' , B q and A q are shown separately. 

lOGeV 2 , namely by 15% to 20% for |t|>20GeV 2 . Comparing with the good fit we obtained with 
(|2"3*)) one might conclude that the data prefers a behavior f q (x) ~ (1 — x) over f q (x) ~ (1 — x) 2 at 
x — > 1, but this would be mistaken as we shall see below. 

In search of a more adequate profile function we demand that 

• the low-x behavior of f q (x) should match the form ()13|) . where we now impose the value a' = 
0.9 GeV 2 from Regge phenomenology, 

• the high-x behavior should be controlled by the parameter A q in (|15|) and not by a' , 

• the intermediate x-region should smoothly interpolate between the two limits, with a few addi- 
tional parameters providing enough flexibility to enable a good fit to the form factor data. 

We found these requirements to be satisfied by the forms 

f q {x) = a'(l - x) 2 log - + 5,(1 - xf + A q x(l - x) (28) 

and 

f q (x) = a'{l - xf log - + B q (l - x) 3 + A q x(l - x) 2 , (29) 

which respectively correspond to n = 1 and n = 2 in (|15|). At large x, the individual terms behave 
like a'(l — x) n+2 , B q {l — x) n+l and A q (l — x) n , which in particular prevents the term with a' from 
being too important in the high-x region. 

As we see in Tables and |H1 a fit to either (|28|) or (|29|) with A u , Ad and B u = Bj as free parameters 
provides a good description of the form factor data. We will comment on setting B u = Bj in Sect. 13.41 
In Fig. 2] we show the profile functions divided by log(l/x) obtained in the fit to ()29|) . as well as the 
individual contributions from the terms with a' , B q and A q . Our criterion that the profile function 
should be controlled by a' at low x and by A q at high x is now well satisfied. To a lesser degree this 
also holds for the fit to (|28|) . where the contribution of the a' term to f u (x) is about 30% at x = 0.7 
and 15% at x = 0.8. The sensitive region of x in both fits is essentially the same as for the fit to 
(|23[) . with the values of x max (t) and (x)t differing by less than 2% from those shown in Fig. [2j The 
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Figure 5: The distance d q (x) between struck quark and spectators, evaluated for the best fits to (|2*%|) 
(dashed line) and (solid line) with = 2GeV. The smallest x value plotted is 5 • 10~ 3 . Shaded 
bands indicate the 1 a uncertainties of the fits as explained in App. [Bj 



difference of x m m(t) in the different fits is more pronounced at small t but below 5% for —t above 
lGeV 2 . 

We see from Tables [7| and |S1 that fd(x) has clearly larger errors than f u (x). This is a feature of all 
our fits (except when we force f u (x) and fd(x) to be equal) and can readily be understood from the 
sum rules (j2J). Due to the charge factors, the integrand giving Ff is dominated by the contribution 
from u quarks. This trend is enhanced by the fact that u quarks are more abundant in the proton 
than d quarks. For large x the ratio d v /u v of parton densities becomes very small indeed (see Fig. EJ), 
and one can expect this trend to persist for H^/Hy at least over some range in t. The combination 
of Ff and Ff provides sensitivity to d quarks, but data on both form factors is only available in a 
relatively small interval of \t\. Improved data on Ff in a wider range of t would be highly welcome 
in this context. 

In Fig. Owe show the distance d q (x) obtained with our fits to (|28|) and to H29|) . For u quarks the 
results of the two fits are fully compatible within their errors up to about x = 0.75. The x-region 
where they differ significantly is outside the range where our fit to the form factors can constrain 
them. For d quarks the results start to differ at lower values of x, but their errors are significantly 
bigger as well. We conclude that the limiting behavior of d q (x) for x — > 1 cannot be determined by 
data on Ff up to |i| around 30GeV 2 . We note that the description of Ff at high |t| is slightly better 
for the fit to (|29j) than for the one to ()28|) . where Ff (t) falls off a bit too fast. One should however 
not interpret this as a preference of the data for n = 2 rather than n = 1 in the power-law falloff (|15|) 
since the situation is opposite for the fits to IJ2HJ) and (|27|). see Tables |H1 El an d El Which value of n a 
fit prefers thus depends on the remaining functional dependence of d q (x). Without data constraining 
d q (x) for x above 0.8 we cannot determine its behavior around x = 1. 

We also see in Fig.[S]that d u from the fit with n = 1 takes values one may suspect to be unphysically 
large only for x above 0.9 or so, where even the forward parton densities are barely known. Similarly, 
for models obtained with the Gaussian wave functions (|14j) and the parameters used in |391 [3] , the 
distance d q stays below 1 fm for a;<0.9. In the kinematic range where these models have been used 
to describe or predict data we hence do not find them physically inconsistent. In the following we 
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will however take the form with n = 2 in Q15j). since its limiting behavior for x — > 1 is more plausible 
than the one with n = 1. An exponent n above 2, which results in a vanishing d g for x — > 1, may 
also be possible |37j . Since the form factor data cannot determine n we refrain however from further 
investigation of this point. We henceforth refer to the fit to Q11JI and (|29|) at fx = 2 GeV as our "default 
fit". Its parameters are 

A u = (1.22 ±0.02) GeV" 2 , A d = (2.59 ± 0.29) GeV" 2 , 

B u = B d = (0.59 ±0.03) GeV" 2 , (30) 

and a' = 0.9GeV~ 2 , with full details given in App. El Note that in this fit the parameter A q has a 
simple physical interpretation as the limit of d 2 /4 for x — ► 1. To a good approximation it also gives 
the value of this quantity over a finite range of large x. 

3.3 Features of the default fit 

In Fig. El we compare the form factor data with the result of our default fit, whose parametric 
uncertainties are shown by the shaded bands. To have a clear representation of the data at large t 
we have scaled the form factors with t 2 . The quality of the fit at smaller t is better seen from the 
"pull" , defined as [iq (data) /F 1 (fit) - 1] and shown in Fig. □ Our fitted GPDs describe Ff within 
5% for —t up to 27 GeV 2 , only the last data point at —t = 31.2 GeV 2 has a larger pull. Apart from 
the data point at —t = 0.255 GeV 2 with its huge errors, our fit describes F™ within 25%. A detailed 
inspection reveals that a large contribution to % 2 is due to five data points in the sample of |U, with 
0.18 GeV 2 < -t < 0.86 GeV 2 . The relative errors for these data are only about 1%, which is below 
the accuracy we are aiming at. We are hence not worried by the comparatively high \ 2 value found 
in our fit. 

The distance d q (x) between the struck quark and the spectators shown in Fig. Elis one of our main 
physics results. Our analysis provides the first data-driven determination of this important quantity. 
It exhibits a significant decrease when going from small to intermediate x. We will comment on the 
increase of dd(x) for x > 0.4 at the end of Sect. 13.41 

In Fig. [31 we show the square root of the average impact parameters (b 2 ) x for the sum and the 
difference of u and d quark distributions, defined as in ® with HI replaced by il" ± H%- That 
(b 2 ) x comes out somewhat bigger for u ± d than for u — d reflects that (b 2 ) x < {b 2 ) x in our fit. The 
points shown in the figure are results from an evaluation of moments j\ dx x m [H u (x, t) ± H d (x, t)] in 
lattice QCD |llj . A quantitative comparison of the two results must be made with due caution. For 
one thing, the values of x in the lattice calculation have been estimated from the ratios of successive 
moments in x at t = 0. We could of course avoid this problem by directly comparing our results for x 
moments with those obtained on the lattice. More importantly, however, the lattice calculation was 
performed for a pion mass of 870 MeV, and an extrapolation to the physical pion mass has not been 
attempted in ^Tj . Indeed, the falloff in t of (F" ± Ff) and (F" — Ff) obtained in that calculation 
|10l 111] is too slow to correctly describe the data for Ff and F™ . Together with the uncertainties 
inherent in our phenomenological extraction, we nevertheless find the overall agreement of the results 
in Fig. |H1 remarkable. 

In Fig. |H1 we also show the profile functions f u and fd themselves. The scale controlling the 
decrease of H%(x, t) with t is seen to depend very strongly on x, not only in the regions of very large 
or very small x but also when going from, say, x = 0.1 to x = 0.3. 

3.4 Variations of the fit 

The 1 a errors quoted in our tables and the associated error bands in the figures reflect the uncertainty 
on the free parameters in a given parameterization of the GPDs, but not the uncertainty due to the 
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Figure 6: Results for the Dirac form factor of proton and neutron with our default fit, defined by 
pip and (|29p with the parameters ()3(jp and the valence quark densities evaluated at scale \x = 2 GeV. 
The solid lines represent the best fit, and the error bands represent the 1 a uncertainties of the fit as 
explained in App. [B] Details on the form factor data and their errors are given in App. EJ 



pull if 




5 10 15 20 25 30 35 



-t [GeV 2 



0.50 

0.25 

0.00 
-0.25 
-0.50 ■ 
-0.75 ■ 



-1.00 



1 



pull F ] 



n 



-t [GeV 2 



Figure 7: The pull [F± (data) fF\ (fit) — 1] for our default fit to the Dirac form factors of the nucleon. 
The error bands represent the 1 a fit uncertainties relative to the best fit. 
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Figure 8: Left: Square root of the average impact parameters (b 2 ) x for the sum and the difference of 
u and d quark distributions obtained in our default fit, compared with lattice QCD results from jllj . 
The smallest x- value plotted is 5- 10~ 3 . Right: The profile functions /„ and fd obtained in our default 
fit (|29|) (solid) with fi = 2GeV. For better visibility fd has been scaled by a factor 1/2. The smallest 
x-value plotted is 2.5 • 10~ 2 . 



choice of parameterization. In order to obtain a better feeling how significant various features of our 
default fit are, we have tested their stability under several modifications of the fit. 

To investigate the difference between the profile functions for u and d quarks we have fitted the 
data to the form Q29|) with or without setting B u = Bd and A u = Ad, see Table |S1 The fit with all 
four parameters free leads only to a slightly better description of the data than our default fit with 
B u = Bd- The fitted parameters are compatible with those in the default fit at 2 a level, but the 
errors on Ad and Bd are much larger. We find that the presently available data on do not warrant 
two free parameters for d quarks. Note that although B u > Bd in this fit, the resulting function 
f u ( x ) only becomes larger than fd(x) for x < 0.1, where the difference between the two functions is 
at most 5%. Already at x = 0.3 the ratio fd/ fu has grown to a value of 1.4, to be compared with 1.3 
in our default fit. A fit where we impose both A u = Ad and B u = Bd cannot adequately describe the 
neutron data. The x 2 i n this subsample is 74 for 8 data points, and the fit result for \F™\ undershoots 
the data by at least 30% for -t > 1 GeV 2 . The same happens if we take the profile function (|28[) with 
fu( x ) = fd(x), see Table A three-parameter fit to (|29|) with the constraint A u = Ad finds B u < Bd- 
It still undershoots the data on \F™\ for —t > lGeV 2 , although not as badly as the fit where both 
A u = A d and B u = B d . 

We conclude that if we insist on having a good description of both proton and neutron form 
factors, we need f d (x) > f u (x) at moderate to large values of x. This implies that the suppression of 
d compared with u quarks seen in the forward parton densities at high x becomes even stronger for 
H% and fljf as \t\ increases. Note that the observed rise of the form factor ratio —F™(t)/Ff(t) implies 
that the flavor contribution Ff(t) decreases faster with \t\ than Fi(t). This is seen by writing 

2F?(t)_ l-n(t) 2F?(t) 

Rl(t) --im-^IW)' l( ) " ■ (31) 
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Figure 9: The ratio d v (x)/u v (x) obtained with different sets of CTEQ6M parton distributions |2()j 
at scale fi = 2 GeV. The largest plotted value of x is 0.97. 

With the data giving R\ ~ 0.37 at — t = 1 GeV 2 one has r\ ~ 0.7, which is clearly different from 
the value r\ = 1 at t = 0. To have r\(t) decreasing sufficiently fast, the damping factor f q in the t 
dependence must be bigger for H% than for iT" when we take the exponential form The same 

trend is observed with the power-law dependence on t we investigate in Sect. 13.61 see Table ITT1 

It is well known that as x becomes larger, the parton densities extracted from data become 
more and more uncertain. In their analysis [20] CTEQ provide 40 sets of parton densities which 
reflect variations from their best fit results that are allowed within errors. Among these we find 
sets 17, 18 and 35, 36 to provide the largest deviations from the best fit parton densities at large x, 
reaching 20% for u v at x = 0.9 and for d v at x = 0.7 and growing well beyond 100% for d quarks 
at higher x. In Fig. 03 we show the corresponding ratio d v (x)/u v (x), which is especially important 
for the simultaneous description of the form factors Ff and in our fit. Repeating our default fit 
with these error distributions as input, we find good stability of the obtained GPDs, see Table 03 
The profile functions f u for the error distributions deviate by at most 3% from the one for the CTEQ 
best fit, as well as fd for sets 17 and 18. With both sets 35 and 36 the ratio of fd obtained for 
the error distribution and for the CTEQ best fit grows from 1 to about 1.1 as x rises from to 1. 
The uncertainties on the forward parton densities thus hardly affect our extraction of the impact 
parameter profile of parton distributions as a function of x. 

We have finally allowed the value of a' in (|29[) to be selected by the fit, and in addition we 
have taken the forward distributions in our ansatz (jll|) at different scales fj,. The results are given 
in Table 1101 In all cases we obtain a rather good description of the form factor data, although 
there is a tendency for the fits to become worse for larger /i. We see that our parameterization of 
GPDs is reasonably flexible to cover a range of factorization scales. This also validates the analytical 
considerations in Sect. 13. 1| which showed that in selected regions of x and t our functional form of 
GPDs is approximately stable under a change of [i. The profile functions f u and fd obtained in our fits 
decrease with /i, precisely as we expect from the evolution equation (jHJ) for (b 2 ) x . In Fig. E3 w e show 
the change of parameters a', A u and B u = Bd with [i. The uncertainty on Ad is too large to observe 
a clear evolution effect. The mild fi dependence in the central values of a' does not contradict our 
general analysis in Sect. 13.11 where a' describes the behavior of the profile function at very small x, 
whereas the parameter a' in our ansatz for f q (x) is relevant in a finite x interval (see Fig. Ej. Given 
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Figure 10: Parameters and their errors obtained in a fit to (jllj) and (|29j) with -B u = i?^ and forward 
parton distributions taken at different scales /i. The corresponding parameters are given in Table fTUl 

this and the uncertainties on a' within Regge phenomenology, we conclude that the fitted values of 
this parameter confirm our assumption that the valence GPDs at small x and t can be described 
by the leading Regge trajectories known from the phenomenology of soft hadronic interactions. The 
results of this exercise also justifies the choice of fixing a' in our default fit, although it is clear that 
its parametric errors underestimate the uncertainty of the profile function in the small- x limit. 

To illustrate the uncertainties of our fit results due to the choice of parameterization, we show 
in Fig. ^2 the average distances d u {x) and dd(x) obtained with selected fits, which provide good 
descriptions for both proton and neutron data. We see that for u quarks there is little influence of 
the parameterization up to x ~ 0.8. For d quarks the uncertainties are larger, due to the lack of good 
neutron data at higher values of t. The curve with the lowest values of dd{x) in the figure belongs to 
the fit in Table |H] where A u = but B u ^ B^, which provides only a moderately good description 
of the neutron data. This is the only fit among those shown where d u and dj differ by less than 10% 
over the entire x range. In all other cases we observe in particular that dd(x) rises for x above a 
certain moderate value, in order to accommodate a rise of the ratio fd/fu- 

3.5 Large t and the Feynman mechanism 

In this subsection we will see that our parameterization of GPDs enables the Feynman mechanism at 
large t, where the struck quark carries most of the nucleon momentum and thus avoids large internal 
virtualities of order t. Let us first consider H->(x,t) in the limit of large x and take a closer look 
at the scaling of momenta, which we denote as shown in Fig. 112b . Defining light-cone coordinates, 
= (v° ±u 3 ) / '"v/2 and vj_ = (v , v 2 ) for any four-vector v, we have x = k + /p + . We choose a reference 
frame where A + = (i.e. £ = 0) and p_|_ = 0. Then t = — , and the virtuality of the active quark 
before it is struck can be written as 

I 2 -4- 1 2 

k 2 = xm 2 + l 2 - 1 —^, (32) 
1 — x 

where m is the nucleon mass. For small (1 — x) we distinguish two momentum regions according to 
the virtuality and transverse momentum of the spectator system: 

soft region: \l 2 \ ~ l\ ~ A 2 , |A; 2 | ~ A 2 /(l — 

ultrasoft region: \l 2 \ ~ l\ ~ (1 - x)A 2 , |/c 2 |~A 2 , (33) 
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Figure 11: The average distance d q (x) for u and d quarks for selected fits. The solid lines show our 
default fit, and dashed lines the fits in the second and third rows of Table |H] and in the second row 
of Table 1101 Dotted lines show the fits in the first row of Table |H1 and the first row of Table Q where 
d q (x) ~ (1 - x)- 1 ' 2 for x-*l. 



where A is a typical scale of soft interactions. Our nomenclature is similar to the one in recent work 
on soft-collinear effective theory |42j . Note that in the "soft region" the spectator partons are soft, 
but the struck quark is far off-shell (in the parlance of [Hj it would be identified as "hard-collinear" ) . 
In the "ultrasoft region" the spectator system has virtuality and squared transverse momentum much 
smaller than A 2 . Such momentum regions do appear in the perturbative analysis of graphs if quarks 
are treated as massless (see e.g. |4.'-}| I44j ) . but one may suspect that due to confinement they cannot 
be important in physical matrix elements. 

An analogous classification holds with respect to the virtuality of the active quark after it is struck, 
with 

I 2 + (l± - (1 - x) A ± ) 2 



k 



xm + V 



1 



x 



(34) 



Note that with our choice of frame, l± — (1 — x)Aj_ is the intrinsic transverse momentum of the 
spectator system in the outgoing nucleon (see e.g. [5]). 

That the struck quark in the soft region is far off-shell is the basis of a perturbative analysis, 
which has long ago been given for the x — > 1 limit of parton distributions or deeply inelastic structure 



(a) 



(b) 



(c) 



k , . k'=k+A 





Figure 12: (a) Momenta of nucleon, active quark and spectator system in a GPD. (b) Perturbative 
mechanism for large x. (c) Hard-scattering mechanism |47j for large t. 
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functions (see e.g. j^Sl) and recently also for the x — > 1 limit of GPDs at arbitrary fixed t [101 • ^ * s 
based on graphs as the one shown in Fig. 112b . where a configuration of three quarks with momentum 
fractions of order | and virtualities of order A 2 turns into a configuration with a fast off-shell quark 
and two soft quarks by successive emission of gluons, which need to be off-shell too. Standard 
perturbative power counting for these graphs in the momentum region just stated gives a behavior 
H$(x, t) ~ (1 — x) 3 for x — ► 1. Their actual calculation in perturbation theory leads however to severe 
divergences dZ 2 /Z 4 in the infrared unless the quark mass is kept finite, which indicates that standard 
hard-scattering factorization [37j does not provide an adequate separation of short-and long-distance 
physics for the mechanism represented by the graphs. While the general power-counting argument 
might still give the correct answer, further details such as the overall normalization are currently 
beyond theoretical control. Resummation of radiative corrections into Sudakov form factors can 
give a stronger suppression than the power-law (1 — x) 3 obtained from fixed-order graphs, but given 
the above difficulties the detailed form of these corrections (let alone their quantitative impact) is 
unknown. We recall that the single logarithms resummed by DGLAP evolution also modify the power 
of (1 — x), see e.g. |ig] . 

Phenomenologically, the limiting behavior of parton densities for x — > 1 is poorly known. Their 
extraction from hard processes is increasingly difficult in this limit due to higher- twist contributions, 
and leading-twist analyses can use data only for values of x where such contributions are under 
control. It is then difficult to infer a power behavior for x — > 1, as our attempt to extract the large- x 
asymptotics of the profile function in Sect. 13.21 has taught us. The powers (1 — x)@ appearing in 
many parameterizations of parton densities are to be seen as parameters describing the approximate 
behavior of these functions over a certain range of large x. The CTEQ6M distributions we use in our 
analysis have powers (3 U ~ 2.9 for u v and (3d ~ 5.0 for d v in the parameterization at the starting scale 
[j, = 1.3 GeV. We find that at fi = 2GeV these distributions are described by u v (x) ~ (1 — x) 3,4 and 
d v (x) ~ (1 — x) 5,0 within 5% for 0.5 < x < 0.9. Taking different x-intervals, these parameters slightly 
change. 

In form factors at large t the relevant x-values of the corresponding GPDs are selected by the 
dynamics. Demanding k 2 and k' 2 to be of the same order, the soft and ultrasoft regions are identified 

as 

soft region: 1 — x ~ A /\/\t\, \k 2 \, |/« /2 | ~ Ai/jij, Z + , Z~, Zj_ ~ A, 

ultrasoft region: l-x~A 2 /|t|, \k 2 \, \k' 2 \ ~ A 2 , Z+, l~ , lj_ ~ A 2 /y/W\- ( 35 ) 

The contribution from the soft region has been analyzed in perturbation theory in the same way as for 
parton distributions |47| I45j . Power counting for the graph in Fig. I12b gives Fi(t) ~ t~ 2 . This is the 
same power as obtained with the standard hard-scattering mechanism shown in Fig. ll2b . where parton 
virtualities are of order \t\. Again one may expect a further damping in the soft region by Sudakov 
corrections, leaving the hard-scattering mechanism as the dominant contribution at asymptotically 
large —t. 

Let us now investigate the large-f behavior of the form factors obtained with our parameterizations 
of GPDs. We see in Fig.|2]that at large t the dominant contribution to comes from a rather narrow 
region of large x. To simplify the analysis we take the large-x approximations q v (x) ~ (1 — x)@ q 
and fq(x) ~ A q (l — x) n . At sufficiently large t the integral Ff(t) = f dx H^(x,t) can then be 
evaluated in the saddle point approximation, obtained by minimizing the exponent in H$(x,t) = 
exp[(3 q log(l — x) + tf q (x)] with respect to x. We then find 

Ff (t) ~ |t|-a+A>/» , l-x.= l^A q \t\) , (36) 
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where x s is the position of the saddle point. We see that for our default fit with n = 2 the dominant 
values of x in the form factor are in the soft region. As observed in |49| . the power behavior in (|36[) 
with n = 2 corresponds to the Drell-Yan relation |5()| between the large-i behavior of form factors 
and the large- x behavior of deep inelastic structure functions. Indeed, the kinematical assumptions 
made in the work of Drell and Yan correspond to the dominance of the soft region. Putting (3 q = 3 
as obtained from dimensional counting we recover the t~ 2 behavior mentioned above. With the 
phenomenological values (3 U = 3.4 and (3d = 5.0 for our input parton distributions at large x one finds 
that for large t the form factor Ff obtained in our fit should fall slightly faster than t -2 , whereas 
Ff should decrease much more strongly. At large t both proton and neutron form factor should then 
be dominated by Ff. Our fit result for Ff in Fig. |H1 shows the large-i behavior expected from these 
arguments, which means that the above approximations are indeed applicable in kinematics where 
there is data. We will show Ff and Ff separately in Sect. 16.21 

Taking n = 1 for the large- x behavior of the profile function, the dominant values of x in the 
saddle point approximation are from the ultrasoft region and give 

F i(t) ~ \t\~ 1 ~ (3q . Our fits to j2H|) 

and Q28JI give a good description of the data for Ff, which are clearly incompatible with such a strong 
decrease in t. This means that the asymptotic behavior has not yet set in at —t ~ 30GeV 2 for these 
parameterizations of GPDs. To understand this, we notice that the approximation f q (x) ~ A q (l — x) n 
works rather well in the interval 0.6 < x < 0.9 for our default fit, but not so well for our fits with 
n = 1. Any inaccuracy in f q (x) appears however exponentiated in H-> and thus in the form factors. 
The validity of asymptotic expansions like (|36j) must hence be carefully investigated on a case-by-case 
basis. 

To quantify how close our parameterizations are from the scaling laws in ()35|) . we start with (x)t 
defined in (|25fl and introduce the quantity 

6 eB (t) = t± log[l-(x) t ], (37) 

which is ^ in the soft and 1 in the ultrasoft region. Using the saddle point approximation for both 
/ dx H^{x,t) and / dx (1 — x)H$(x,t), one readily finds that 8 e g(t) tends to n~ 1 at asymptotic values 
of t. In Fig. El we show 5 e s(t) for our default fit and find that for —t above 5 to lOGeV 2 the scaling 
of the relevant x values in the form factor integral is indeed the one for the soft region. In the same 
figure we show 

Aeff(t) = [1 - (x) t ] yj\t\, (38) 

which according to (|35|) should be of typical hadronic size for t where the soft scaling law applies. 
This is indeed the case for our fit. We remark that with the saddle point approximation we find 
A c g = 1.26 GeV when neglecting the d quark contribution to Ff at large t and taking the parameters 
A u = 1.22 GeV -2 and (3 U = 3.4 for u quarks. This shows once more the relevance of asymptotic 
considerations for our default fit in kinematics where data is available. 

The values of 5 e R for our fit to (|28|) . where n = 1, differ from those shown in Fig. El by at most 
8% and the values of A e g- by at most 4%. In the large- i region of present data the soft contribution 
hence also dominates for this fit, where we find that ultrasoft behavior with 5 e fi = 1 only sets in 
for \t\ well above 100 GeV 2 . We remark that in our previous work [5] we used GPDs obtained from 
Gaussian wave functions (|14|). which had a profile function decreasing like (1 — x) at large x and 
gave an asymptotic behavior F\(t) ~ t~ 4 . The power counting for the "soft overlap mechanism" we 
set up in that work corresponds to the ultrasoft region in the parlance of our present paper. The 
considerations of this section make it clear that this (physically suspect) region is not the dominant 
one in the kinematics where the model of [2] has been used for phenomenology. 

We conclude that the description of Ff provided by our fitted GPDs supports the hypothesis 
that for —t from about 5 GeV 2 to several 10 GeV 2 the dynamics is dominated by the Feynman 
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Figure 13: Dynamical interpretation of our default fit. Left: The effective power S e g, which describes 
the scaling of (1 — x) as a function of \t\. Right: The effective soft scale A e g- as a function of \t\. 



mechanism in the soft kinematics specified above. We emphasize that by itself our result cannot 
exclude the dominance of the standard hard-scattering mechanism at large t. Applied to GPDs |5J 
this mechanism results in a behavior 

H%(x, t) ~ t~ 2 h(x) for -t ->• oo at fixed x (39) 

up to logarithms in t. The function h(x) diverges for x — > 1, signaling the breakdown of the fac- 
torization scheme in that limit, see |52j . Our parameterization (jllj) does not tend to the factorized 
form (|39j) in the large-i limit and thus does not incorporate the physics of the hard-scattering mech- 
anism. The dominance of this mechanism in Ff at experimentally accessible t is very doubtful: to be 
close to the data one requires proton distribution amplitudes for which a substantial fraction of the 
form factor comes from configurations where partons are soft and the approximations of leading-twist 
hard-scattering factorization are inadequate. References can e.g. be found in Sect. 10 of I n ° ur 
present work we make the assumption that in the t region we consider, the hard-scattering mechanism 
is not dominant and that the Feynman mechanism controls form factors at large t, despite its possible 
Sudakov suppression in the asymptotic limit. 



3.6 The t dependence 

In this subsection we explore an ansatz for the t dependence of H% that is different from the exponential 
form we have used so far. To cover a range of possibilities we take 

H*(x,t) = q v (x)(l- t -fjj^J \ (40) 

with different powers p. In the limit p — > oo we recover the exponential (|11|) . The corresponding 
impact parameter distribution q v (x,b 2 ) can be expressed in terms of the modified Bessel function 
AT p _i and satisfies positivity. At fixed x the form (|4T))) gives a power-law falloff H$(x,t) ~ \t\~ p for 
\t\ — > oo. 
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Before proceeding to fits let us investigate some general properties of this ansatz. For p ^ 1 the 
form (jiU)) is not stable under DGLAP evolution. To see this we observe that it satisfies 

(|^^(M))(^log23(M)) =1 (41) 

for all x and t, and starting from (J7J) calculate the evolution of the left-hand-side of this relation. If 
at a given scale [i the GPD satisfies (|41jl then this evolution equation simplifies to 



.2 A 



» ,.,2 



X 



where for better legibility we have omitted the arguments t and fx in the GPDs. Except for trivial 
functions H%(x) the right-hand-side of this expression is positive and furthermore depends on x, so 
that relation ()4*T)) is destroyed by evolution. 

At very large t the form factors obtained with (|4Uj) are again dominated by large x provided that 
pn > (3 q , where as before we assume the large- x behaviors q v (x) ~ (1 — x)@ q and f q (x) ~ ^4 g (l — x) n . 
For u quarks and n = 2 this condition is fulfilled even when p = 2. One can then again use the saddle 
point approximation and finds 

Ff(t) ~ \t\-^MI-, X -x s -)A q \t\) ' , (43) 



where remarkably the large-t behavior of F± is independent of p. The dominant x in the form factor 
integral are in the region of the soft Feynman mechanism discussed in the previous subsection. Note 
that in this region tf q (x) is of order 1, so that the approximation giving a power-law H$(x,t) ~ \t\~ p 
is not valid. 

We note that the form 1)40(1 does not have the Regge behavior 1)10)1 at small x and t. Depending 
on p, this can be significant in kinematics appropriate for Regge phenomenology. Taking for example 
x = 10" 3 and -t = 0.4 GeV 2 , we have tf(x) « a't log(l/x) « -2.5 with a' = 0.9 GeV" 2 , and 
(1 + 2.5/p)~ p is about twice as large as exp(— 2.5) if p = 3. Having lost the connection with Regge 
phenomenology, we will not impose a fixed value of a' in our fits. A logarithmic increase of (b 2 ) x at 
small x seems however plausible from a more general point of view, and we keep the analytic form 
of our profile function as in our default fit. One could of course construct parameterizations which 
interpolate between an exponential t dependence at very small x and a power-law when x becomes 
larger. Since the small-x region does however not affect our fit to nucleon form factors for —t above 
1 GeV 2 or so, we shall not pursue this possibility here. 

The results of our fits to ()40)) and ()29|) are shown in Table for selected values of p, where 
for easier comparison we have also given the result of our exponential fit with free a' discussed in 
Sect. 13.41 As p decreases from oo, the quality of the fit first improves (except for the neutron data). 
The smallest x 2 is attained for p = 2.5, and further decrease of p makes the fit worse. With p = 2 no 
good description of the large-i data can be achieved: the fit result for Ff systematically overshoots 
the data above -t = 9 GeV 2 , by about 15% when -t > 19 GeV 2 . We find the same qualitative 
picture when fixing a' = 0.9 GeV -2 in the fit. The lowest % 2 is then obtained at p = 3, and for p = 2 
the description of the large-t data is again very bad. It is important to realize that, despite a clear 
decrease in \ 2 ■> the description of the data for Ff does not improve dramatically when going from 
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Figure 14: Left: The distance d u (x) between struck quark and spectators obtained with fits to the 
power-law (jjOJ) with p = 2.5 and p = 3, compared with the result of our default fit, which corresponds 
to p — > oo and differs only slightly from the corresponding fit with free a'. Right: The results for 
H%(x,t) at -t = lOGeV 2 from the same fits. 

p = oo to p = 2.5. Except for the two data points with the largest values of —t and the data point 
with —t = 0.86 GeV 2 (which has the highest —t in the sample of j41j). the pull for Ff of our fit with 
p = 2.5 is at most 3% in magnitude, whereas it is at most 3.5% with our default fit. 

With decreasing p the profile functions f q {x) increase significantly, which is shown in Fig. 1141 for 
u quarks. This can be readily understood: the GPDs in (1401) decrease much smaller in the variable 
\t\f q (x) for small p than for large p, so that small p requires a larger f q (x) in order to not overshoot 
the form factors at large \t\. 

In the figure we also show the x dependence of H% at large t. Its shape is much broader for the 
fits with smaller p and reaches down to smaller x. This reflects that when p becomes smaller it takes 
significantly larger t to suppress H$ at a given x. For small p the asymptotic behavior Q43JI obtained 
with the saddle-point approximation should hence set in at much larger t. Taking — t = 40 GeV 2 we 
find indeed that (x)t is only 0.6 for p = 3 and 0.5 for p = 2.5, and that 6 e s only becomes 0.3 for p = 3 
and 0.15 for p = 2.5, to be compared with the corresponding values (x)t = 0.8 and S e g = 0.5 for the 
default fit. 

As is seen in Fig. 1141 the GPDs with power-law form (|4U|) do not vanish for x — > at large t, 
so that small x also contribute in the high-t form factors to some extent. This is in contrast to the 
GPDs with exponential t dependence. Using that q v (x) ~ x~ a for x — ► 0, one readily finds that 
H$(x,t) ~ x~ a ~ a t vanishes in this limit as soon as —t>a/a', i.e., already for —t above 0.7 GeV 2 
or so. 

The fits we have shown make it clear that a significant ambiguity remains when one tries to 
determine the correlated dependence on x and t of GPDs from experimental knowledge of their 
integrals over x and of their limits at t = 0. Without data on observables depending on both 
variables we need a certain theoretical bias in our extraction of GPDs. We regard the connection 
of our default fit with Regge phenomenology at small x and t and with the dynamics of the soft 
Feynman mechanism at large t as physically attractive features and retain the exponential form pijl 
for the t dependence of H%(x,t). 
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4 The polarized distribution H v (x,t) 

In this section we investigate the distribution of polarized quarks and its connection with the isovector 
axial form factor Fa(£) of the nucleon. The relevant sum rule reads 



F A (t)= J 1 dx (H u (x,t,fx 2 )-H d (x,t u 



/OJ. (44) 

Its value at t = is given by the axial charge Fa(0) and well known from (5 decay experiments. 
The measurement of Fa{£) covering the largest t range [53] nas been performed in charged current 
scattering n — > fi~p. The result has been given in the form of a dipole parameterization 



F A (0) 
{l-t/Mlf 



with F A (0) = 1.23 ± 0.01 and M A = (1.05 tg^) GeV for a measured range 0.1 GeV 2 < -t < 3GeV 2 . 
The resulting t dependence of FA(t) is not too different from the one of the proton Dirac form factor: 
we find that Ff(t) can be approximated to 9% accuracy for —t < 3 GeV 2 by (1 — t/Mp)~ 2 with 
Mc = 0.98 GeV. Note that the well-known dipole mass of 0.84 GeV does not refer to Ff but to a 
dipole parameterization of the magnetic form factors of proton and neutron, G P M and G 1 ^. 

Many other measurements of FA(t) in either charged current scattering or ep — > e7r + n have also 
been presented in terms of a dipole mass, with a considerable spread of results for — t < 1 GeV 2 (see 
e.g. |54|). We remark that while a dipole form (|45|) provides a simple and compact parameterization, 
it is not well suited for a description of data beyond a certain accuracy. It is instructive in this respect 
to perform dipole fits to the data on Ff or on G P M in different ranges of t and to observe the shift in 
the dipole mass. 

Given the data situation for FA(t) we do not attempt a fit of GPDs to the sum rule (|44|). but 
rather test the simple ansatz 

H«(x, t) = Aq v (x) exp [t~f q {x)\ (46) 

with the profile functions f q (x) taken equal to f q (x) obtained in our default fit for the unpolarized 
distributions H%. In physical terms this ansatz assumes that the distribution of quarks minus an- 
tiquarks in the transverse plane does not depend on the quark or antiquark helicity relative to the 
helicity of the proton. For our evaluation we take the polarized parton densities Aq v (x) from |21j . 
more specifically the NLO distributions in their scenario 1 at // = 2 GeV. 

Since the axial form factor has positive charge parity, it is not directly connected with the valence 
quark distributions. Instead we have 



F A (t) = J dx (h%(x, t) - H*(x, tj) + 2 dx (H u (x, t) - H S (x, tj) , (47) 

where the generalized antiquark distributions are given by H^(x,t) = H q (—x,t). The flavor non- 
singlet combination Au(x) — Ad(x) of forward densities is poorly known, and at present there is 
no experimental evidence that it might be large [HS1- To make a motivated ansatz for its analog 
at finite t is beyond the scope of this work. For a very crude estimate we have taken H^(x, t) = 
Aq(x) exp[tf q (x)] with the polarized antiquark distributions from |21j . where Au(x) = Ad(x). The 
resulting contribution from antiquark GPDs in l|47jl is below a percent, both for f q (x) from our default 
fit and from the corresponding four-parameter fit in Table |HJ where B u ^ B^, This is because the 
profile functions for u and d quarks differ mostly at larger x, where antiquarks do not abound. 

The present uncertainties on polarized quark distributions are significantly larger than those on 
their unpolarized counterparts. We have calculated the resulting error on FA(t) using the covariance 
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Figure 15: Axial form factor -Fa(*) obtained from the valence distributions ()46|) with polarized 
densities from at fi = 2 GeV 2 and the result of our default fit for f q (x) specified in Sect. 13.21 The 
contribution from sea quarks has been neglected. The error band gives the 1 a uncertainties of the 
fit to f q (x) added in quadrature to those of Aq(x). The curve and band shown up to —t = 5 GeV 2 
represents the dipole parameterization from with Fa(0) = 1.23±0.01 and Ma = (1.05 iolii) GeV. 

matrix on the parameters in the parton densities provided in [2^. This error is at least a factor 
of 5 larger than the error resulting from the uncertainty in f q (x). For estimating the parametric 
uncertainties of FA(t) obtained with our ansatz, we have added the errors from the two sources in 
quadrature. 

In Fig. we show the contribution to Fa(£) from the valence quark GPDs specified above. 
We compare this with the dipole parameterization of [HS], which we have extrapolated up to —t = 
5 GeV 2 . Our result undershoots the central value of that parameterization by at most 13%, with the 
largest discrepancy for —t between 0.5 and lGeV 2 . It is however consistent within the errors of the 
parameterization in the full range 0.1 GeV 2 < — t < 3 GeV 2 of the measurement in |53| . We note that 
one obtains a larger form factor with the ansatz (|46j) when taking a smaller profile function f q (x). 
This is indeed physically allowed, whereas taking f q (x) > f q (x) would violate the positivity of parton 
densities in impact parameter space (see Sect. I5.2|) . 

In conclusion, we find that the present data on the axial form factor is consistent with valence 
quark dominance in the sum rule (|47jl and with only a weak helicity dependence in the transverse 
distribution of valence quarks at not too large values of x. We note that x max (i) for i^i(t), defined in 
analogy to J2U), equals 0.73 at -t = 3 GeV 2 . 

5 The helicity-flip distribution E from the Pauli form factors 
5.1 General properties 

The GPDs E q are related to the Pauli form factors of proton and neutron through the sum rules 

|^(x,t,At 2 )-^(x,t,At)l , 




if (t) = / dx 
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F-f(i) = f Ax^Eiix^^-lE^x.t,^)], (48) 

<J 

where in analogy to we have introduced valence distributions 

E%(x, t, fj 2 ) = E*(x, t, /i 2 ) + E«(-x, t, ii 2 ) (49) 

for quarks of flavor q in the proton. Contributions from sea quarks cancel in (|48j) . We have neglected 
in (|48|) the contribution from strange quarks, as we did for the Dirac form factors. The scale depen- 
dence of E v is described by the same DGLAP equation as the scale dependence of H v , since both 
distributions belong to the same quark operator. 

The distribution E q (x,t) describes proton helicity flip in a frame where the proton moves fast 
(or more precisely light-cone helicity flip, see e.g. Sect. 3.5 of d]). Since quark helicity is conserved 
by the vector current, proton helicity flip requires orbital angular momentum between the struck 
quark and the spectator system. This becomes for instance manifest by writing E q as the overlap of 
light-cone wave functions, whose orbital angular momentum must differ by exactly one unit. Another 
manifestation is Ji's sum rule for the combination H q (x,t) + E q (x,t), see Section f(i. 31 

E q admits a probability interpretation in impact parameter space if one changes basis from lon- 
gitudinal to transverse polarization states of the proton [21]. More precisely, one considers states of 
definite proton transversity, which is the light-cone analog of transverse polarization, see e.g. |56j . 
The distribution 

b y 

q*{x,b) = q v (x,b) —^e q (x,b) (50) 

m do 

gives the probability to find an unpolarized quark with momentum fraction x and impact parameter 
b = (b x ,b y ) in a proton polarized along the x direction, minus the corresponding probability to find 
an antiquark. Here we have introduced the Fourier transform 

e q v (x,b) = J ^e- ibA E«(x,t = -A 2 ). (51) 

Note that e q (x,b) and q v (x,b) depend on b x and b y only via b 2 . We see from (|5U|) that target 
polarization along the x-axis induces a shift in the quark distribution along the y-axis. As explained 
in this effect is consistent with the classical picture of the polarized proton as a sphere rotating 
about the x-axis and moving in the z-direction (see also (HZ])- The average displacement of this shift 
is 

fd 2 bby q X(x,b) J_ E q (x,0) 
K )x fd 2 bq*(x,b) 2mHS{x,0) 7 1 ' 

and its scale evolution is given by 

V 2 ^l(b V )l = --^T C-Pi-^q^Ub^l-^Yl (53) 



dfi 2 x q v (x) J x z \Z' 

in analogy to the evolution of (b 2 )^. The corresponding shift for the distance between the struck 
quark and the spectator system is 

Sq W = rr|. ( 54 ) 

in analogy to the distance function d q {x) we introduced in @. 

The impact parameter space distributions satisfy inequalities which insure that the quark densities 
for various combinations of proton and quark spins are positive. Using the methods of [SS] one finds 
in particular [1H] 

2 

l {x,b)\ < \q(x,b) + Aq(x,b)]\q(x,b) - Aq(x,b)} , (55) 
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where e q (x,b), q(x,b) and Aq(x,b) are the respective Fourier transforms of E q (x,t), H q (x,t) and 
H q (x,t), defined in analogy to (J3J) and (|51|). The bound <|55|) is stable under leading-order DGLAP 
evolution to higher scales, and a closer look at its derivation shows that it should be valid when \x 
is large enough for the application of leading-twist factorization theorems to the exclusive processes 
where GPDs appear. For a discussion and references see |19j . 

Multiplying (f55|) with b 2 and integrating over b, one obtains after a few steps an inequality for 
GPDs in momentum space |49j . 

(£"(x,0)) 2 < 4m 2 [g(x) + Aq(x)\ \q(x) - Aq(x)] ^ In [H q (x, t) ± H q (x, t)]^ , (56) 

where either the sum or the difference of H q and H q may be taken. The t-derivative on the right- 
hand side is \ times the average squared impact parameter of quarks with positive or negative 
helicity. According to our discussion in Sect. 13. ll this quantity may be expected to decrease at least 
like (1 — x) 2 in the limit x — ► 1. Since in addition the longitudinal polarization |Ag/<7| of quarks is 
phenomenologically seen to grow as x becomes large, the inequality ()56|) severely limits the high-x 
behavior of E q (x,t = 0) and will be an essential input in the following. 

The positivity constraints ()55|) and (|56j) hold for the distribution of quarks (i.e. for x > 0) and 
have analogs for antiquarks. They do not hold for the valence combinations we aim to determine 
in this work, which are differences of quark and antiquark distributions. It is however natural to 
neglect antiquarks at large enough x, which is in fact the region where the bounds give the strongest 
constraints. With this proviso in mind we will in the following use (|55j) and (|56|) for the valence 
distributions. 

Our discussion of the bound (|56j) implies that the average displacement {b y ) q should vanish in the 
limit of x — ► 1. At large enough x it should hence be a decreasing function, which according to (|53j) 
implies that {b y ) q becomes smaller when evolving to higher scales [i. We comment on the small- x 
behavior of this displacement in the next subsection. 

5.2 Ansatz for E v (x,t) 

Our ansatz for E v is taken in analogy to the other GPDs as 

E q {x,t) = e q (x) esxp[tg q {x)], (57) 

where we use the notation 

e q (x) = E q (x,t = 0) (58) 
for the forward limit. The normalization integrals 

/ 

Jo 



dx e q (x, t = 0, /i) = K q (59) 



give the contribution of quark flavor q to the anomalous magnetic moment of the proton (up to quark 
charge factors). Neglecting the contribution from strange quarks one has according to ()48|) 

k u = 2k p + K n w 1.67 , = k p + 2k„ w —2.03 . (60) 

For the profile function in (|57|) we take the same form as in our default fit for H q , 

g q (x) = a'(l-xflog^ + D q (l-xf + C q x(l-x) 2 with a' = 0.9 GeV" 2 . (61) 
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The motivation from Regge phenomenology for the behavior of g q at very small x applies in the same 
way as for f q . Indeed the Regge exchanges contributing to H$ and E% are the same, and only their 
coupling strengths differ. We therefore take the same fixed value of a' as in our default fit for H$. 
The exponential t-dependence for H$ and E$ taken in (|46|) and (|57|) gives 



Aq v (x, b) 



J_ Aq v (x ) 
^ f q (x) 



cxp 



e v (x,b) 



J_ e v (x) 
4vr g q (x) 



exp 



Ag q (x) 



(62) 



in analogy to the form (|12j) of q v (x,b). According to (|55|) one must have |A^(a;, b| < q v (x,b) for all 
6 when antiquarks are negligible, which implies < as anticipated in Sect. With our 

simplified ansatz f q (x) = f q (x), the bounds (f55|) and (|56|) respectively read 
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and 



exp 



(el) 2 < Am 2 f q {q v f - (Aq v ) 



2g q 



< 



(q v ) 2 - (Aq v f 



exp 



Vq 



(63) 



(64) 



when antiquarks can be neglected, where for the sake of legibility we have omitted the argument x in 
all functions. Because of the factor b 2 on the left-hand side of (|63|) we must have 



9q(x) < fg( x ) 



(65) 



with strict inequality, otherwise the bound will be violated at sufficiently large b 2 . Multiplying both 
sides of ()63|) with exp[& 2 /(2/ 9 )] and then maximizing the left-hand side, one finds that the bound is 
most stringent for b 2 = 2g q f q /(f q — g q ), where it reads 



(e2) 2 <4m 2 e 1+k W^ 



(Jq-9q) (<lv?-{Aq v f 



(66) 



We note that the bound (|64|) is weaker than (|66|) but has the practical advantage to be independent 
of the profile function g q . If we require the distance d q (x) between struck quark and spectators in an 
unpolarized proton to stay finite in the limit x — ► 1, then this bound guarantees that the shift s q (x) 



of this distance also remains finite, given that el = 2m (1 
For the shape of e% we make the time-honored ansatz 
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\q v s q and f q = (1 
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e%{x) 



N q K q X~ 



(1 



(67) 



whose analog for q v and Aq v gives a reasonable first approximation of phenomenologically extracted 
parton densities. The factor 

_ T(2-a + P q ) 

- m-aw+M (68) 

ensures the proper normalization ()59|) . In the small-x limit a takes the role of a Regge intercept if 
one assumes dominance of a single Regge pole, see Sect. 13.11 From Regge phenomenology and from 
experience with q v one expects a ~ 0.5 for both u and d quarks. Note however that the form (|67|) 
is sensitive to a over a finite interval of x, and we do not have enough data to introduce further 
parameters which would make the description of el more flexible. For the average displacement (b y )^. 
in 1)52(1 we expect a relatively weak x-dependence at small x, say a power-law x s , where 5 should be 
of order 0.1 but may be positive or negative. In line with our treatment of a' we have taken a single 
parameter a for u and d quarks in all our fits. Trying to determine a flavor dependence at the level 
of what is seen in the distributions q v (x) is beyond the accuracy we can hope for in this study. 
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5.3 Fit to the Pauli form factors 



Before proceeding to the fits of E% we would like to point out two features of the data on the Pauli 
form factors. As we saw in Sect. 13.41 a comparison of the Dirac form factors for proton and neutron 
clearly shows that with growing \t\ the ratio of Ff and i* 1 " must become smaller than its value \ at 
t = 0. Let us investigate the evolution with \t\ for the analogous ratio of individual flavor contributions 
F$(t) = Jg 1 dxE$(x,t) to the Pauli form factor. Writing r 2 = (k^ 1 F^) / (k' 1 Fg) we obtain 

k^F? 1 + 0.71 (r 2 - 1) , 
^ = ^o ^ in ,o , « 1 + 0.33 (r 2 -l), (69) 



K 



v l F% l + 0.38(r 2 -l) 



where we have inserted the values of the anomalous magnetic moments and in the last step expanded 
in the deviation of r 2 from its value at t = 0. In Fig. ^] we show the normalized form factors k" 1 ^ 
and K^Frf, weighted with a factor \t\ to make the region where we have neutron data more visible. 
We observe that the five neutron data points with —t between 0.25 and 1.2 GeV 2 have a tendency to 
be above the data for the proton. We note that these points are from three different measurements 
[591 160j , |61| and , so that this effect is at least not due to a normalization problem in a single 
experiment. According to (|69j) . a positive value of R 2 — 1 implies a positive value of r 2 — 1 about 
three times as large, so that K^ 1 i ? 2 u (t) must decrease faster than n^F^ (t) starting from t = 0. This 
is the opposite of what is found for the Dirac form factors, and it would be interesting to have better 
neutron data to see if this trend is confirmed, and possibly reversed at larger t. With the data and 
errors at our disposal, the trend has a clear effect on our fit of E™ and E%. 

In Fig. El we show -F 2 weighted with a factor t 2 (as we did for the Dirac form factors) and find a 
striking plateau for —t between 2.5 and 5.5 GeV 2 . A behavior F% ~ t~ 2 is certainly not expected in 
the large-t limit. The plot in Fig. El thus instructs us that observables may exhibit an approximate 
power-law behavior in an intermediate range of a kinematical variable, which has little to do with the 
asymptotic behavior. 

Apart from their connection with the Pauli form factors, we have at present no phenomenological 
constraints on the shape of the GPDs E q (x, £, t). With data for F 2 going up to —t = 5.5 GeV 2 for the 
proton and even less for the neutron, we must expect that a significant range of functions e%{x) and 
g q {x) in (|57|) is able to describe the form factors. In particular, the fit can partially accommodate a 
decreased value of (3 q by simultaneously increasing g q , thus partially compensating a shift of e% (x) to 
larger values of x by a stronger suppression through the exponential factor exp[tg q (x)] in that region 
of x. The bounds (fHljl and provide lower limits on j3 q and upper limits on g q {x). In the opposite 
direction we have only the requirement that g q {x) must be positive for the exponential ansatz (|57[) 
to make sense. 

Performing fits to the data on F 2 and i 7 ^ where the parameters C u and C d are left free, we find 
that the values of Cd and (depending on the details of the fit) also of C u are too large to fulfill the 
bound (|65|). which results in an exponentially strong violation of (|63|) for large enough x and b 2 . We 
have therefore in the minimum-x 2 fits imposed the bounds 

C u < A u = 1.22GeV~ 2 , C d < A d = 2.59 GeV" 2 . (70) 

The values of A u and A d in (|7U|) are those obtained in our default fit to the Dirac form factors (see 
Sect. EH}- Directly implementing (|64|) or ()66|) in the fit would be more involved, and we have instead 
verified these bounds only for the fit results. We remark that with our input parton densities at [i = 
2 GeV (CTEQ6M [20] and the NLO distributions in scenario 1 of [H]) the inequalities |Ag(x)| < q(x) 
and |Ag v (a;)| < q v (x) are violated for the best fit values of the distributions at x above 0.6 or 0.7, but 
are well satisfied within the 1 a error bands on the polarized densities from |22- We therefore require 
the bounds (|63|) . (|64|l. (|66|) to hold within these errors. 
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Figure 16: Data for the normalized form factors K p i 7 ^ and K n weighted with 
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Figure 17: Data for the weighted form factor ftF^ (t). Full symbols: results obtained with the recoil 
polarization method [881 1761 Open symbols: results from the Rosenbluth separation in |9Uj . 
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Notice that when fitting E%(x,t, fi 2 ) to the Pauli form factors alone, we cannot determine the 
scale \x since the form factor integrals in (|48|) are scale independent. We fix [i implicitly in our fit 
by requiring the positivity condition (j66|) to hold for the parton densities at [i = 2 GeV and the 
associated profile function f q (x) obtained in our default fit of H-*. 

To avoid dealing with too many free parameters, we will present fits where we imposed a = 0.55. 
Leaving a free in those fits which give a good description of the data, we obtain values between 0.5 
and 0.6. Following the discussion at the end of the previous subsection, this validates once more our 
assumption that simple Regge phenomenology provides a guide for the behavior of GPDs at small x 
and small t. For those fits which fail to describe the data, making a a free parameter does not help. 
The conclusions we draw in the following are hence not affected by having fixed the parameter a. 

Even then we still have six parameters to determine, namely (3 q , C q and D q for each quark flavor. 
To gain some insight into the range of allowed parameters, we have performed a series of fits with 
fixed values of (3 U and (3d, leaving the other four parameters to be fitted. The resulting values of 
X 2 /d.o.f. are shown in Tableland the results for a subset of these fits in Table IT21 Distributions 
with (3 U = 3 or with (3d = 4 violate the positivity condition (|B^|) at large x and have been discarded. 
A number of observations can be made in these fits: 

1. Both C q and D q increase when (3 U decreases, and for (3 U = 4 the fit selects the maximum allowed 
value C u = A u . For all fits shown in Tabled the fit selects the maximum Cd = Ad, and only 
Dd increases when (3d becomes smaller. 

2. The minimum x 2 is very flat as a function of (3 U and even more of (3d- In particular there is 
only a slight preference to have f3 u < (3d, with good fits being also obtained for (3 U > (3d- 

3. The fitted parameters D u and Dd differ significantly within their errors — much more strongly 
than in our our fits to the Dirac form factors (see Sect. l3~I|) . 

Table 1: Values of x 2 /d.o.f. obtained in fits to (fB7]l. (|6T|) and (JHZI) with fixed a = 0.55. Free 
parameters are D u , Dd and C u , Cd under the constraints (|70|) . 
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These fits with a grid of fixed values (3 U and (3d are not well suited for the propagation of correlated 
errors. We have thus performed a six-parameter fit with free (3 U , (3d, C u , Cd, D u and Dd, fixing only 
a = 0.55. The minimum x 2 is achieved for (3 U = 3.99 and (3d = 5.59, and the parameters C u and 
Cd take their maximum values given in (|70|), To avoid the treatment of errors on parameters at the 
boundary of their allowed range, we have fixed their values to C u = 1.22 GeV -2 and Cd = 2.59 GeV -2 
and repeated the fit. The resulting four-parameter fit still has a very large error of 2.67 on (3d- In 
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Figure 18: Region of x (white region) which accounts for 90% of F^if) in our default fit for E$(x,t). 
The upper and lower shaded x-regions each account for 5% of F^if). The thick line shows the 
average {x)t in the form factor integral. 

addition, this error is strongly correlated with the error on D d , with a linear combination of the two 
parameters being essentially undetermined. The available data on the Pauli form factors do not allow 
us to extract four independent parameters within an acceptable accuracy. As a simple solution one 
might impose the constraint (3 U = f3 d , but this results in a fitted value of f3 q = 4.21, which is in clear 
violation of the positivity bound (|64[). As an alternative we have fixed the difference j3 d — U to the 
value 1.6 obtained in the six- and four-parameter fits just discussed. This value is then to be regarded 
as an external input to a three-parameter fit with 

(3 U = 3.99 ± 0.22, At = A* + 1-60, 

C u = 1.22GeV~ 2 , C d = 2.59GeV" 2 , 

D u = (0.38 ±0.11) GeV" 2 , D d = -(0.75 ± 0.05) GeV -2 , (71) 

and a = 0.55. Full details are given in App. [Bj We take this as our default fit for E% in the rest 
of this paper. We remark that a four-parameter fit where we additionally leave a free finds it to be 
0.55 ± 0.03, with the remaining parameters essentially as in (|71|l. 

The profile function g u obtained with this fit is smaller than f u by at most 10% for all x. In 
contrast, g d is significantly below its counterpart f d at moderate values of x, namely by as much as 
a factor of 2 for x ~ 0.2. This can be understood from the observation we made at the beginning 
of this section. The forward limit e d (x) in this fit is concentrated at smaller values of x than e u (x), 
which favors \F$ | decreasing faster with \t\ than F%- To obtain the reverse trend at small t, which is 
favored by the data, the fit requires a rather weak damping factor g d in the t dependence of E% at 
the relevant x values. In fact, one has g d < g u up to x< 0.45. For larger x the hierarchy is reversed, 
so that eventually \F^\ will decrease faster than i 7 ^ when t becomes large. 

To quantify the sensitivity of our fit to the range of x in the sum rules (|48|) we plot in Fig. 1181 the 
quantities x m i n (t), x max (t) and [x)t for the Pauli form factor of the proton, defined in analogy to (|2*1|) 
and (|25j) by replacing H% with E% and Ff with F%. We find that we are sensitive to values up to about 
x ~ 0.7 with the form factor data at hand. In Fig.^Jwe show the result of our fit in comparison with 
the data, and in Fig. [2Q1 the corresponding pull. With the exception of three (somewhat outlying) 
data points, our fit describes F% within 5% over the entire t range. The description of F£ is of similar 
quality, except for the two data points with the highest —t, where the central value of the fit is just 
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Figure 20: The pull [F2 (data) /F2 (fit) — 1] for our default fit to the Pauli form factors of the nucleon. 
The error bands represent the 1 a fit uncertainty relative to the best fit. 
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compatible with the errors of the data. 
5.4 Flavor structure 

The fit we have just described exhibits a clear difference in the parameters of the u and the d quark 
distributions. To test the significance of this result, we have performed a fit with the same f3 q , C q 
and D q for both flavors, leaving also a as a free parameter in order not to be biased by its particular 
value. The bounds (f7fl|) imply of course C q < 1.22 GeV -2 in this case. Such a fit does not give 
a satisfactory description of the neutron data, with a x 2 of 106 for the 8 data points according to 
Table ITU This fit undershoots all data for \F£\ with —t below 1.5 GeV 2 by 5% to 10%, which is a 
large effects given the errors on the data in that region (see App. ^J. To see whether the positivity 
constraints are responsible for this failure, we have repeated the fit without restricting C q . The fit 
then selects parameters C q and [3 q that badly violate the positivity constraints, but the resulting FQ 
is barely changed for — t < 2 GeV 2 and thus equally inadequate. 

Having seen that the data and errors we use require different k" 1 ^ and K^ 1 E^j, we would like to 
explore for which of our parameters a flavor dependence is required most. From the fits with fixed (3 U 
and Pd presented in Tables ITI and TFZl we have already seen that n~ 1 e^{x) = K^" 1 e^(x) for the forward 
limit is compatible with existing data. Concerning the profile functions, one may ask whether the 
data can be described when taking D u = Dd, as was the case for the analogous fits of H$. A fit with 
free (3 U and (3d and this constraint selects fid much too small to comply with the positivity bound 
(|64[). so that instead we have performed a series of fits with fixed j3 u and (3d, leaving free parameters 
D u = Dd in addition to C u and Cd subject to (|7U|). The resulting values of x 2 /d.o.f. are given in 
Table |2] and details of selected fits in Table 1151 We find that a good description of both proton 
and neutron data is possible only if (3 U is significantly larger than fid- In analogy to the discussion 
at the end of Sect. 15. 3( this can again be understood from the necessity of the fit to accommodate 
K u lp 2 < K d lp 2 at sma11 In Table FTTH we observe that with increasing fi u the fit requires smaller 
and smaller profile functions g u (x), otherwise F% would decrease too fast to be compatible with the 
proton data at larger values of t. For the largest [3 U in the table this does not seem a physically 
very plausible scenario. In the case [3 U = 11 and fid = 5 the fitted values even lead to g u (x) < for 
0.5 < x < 0.8, so that the corresponding entry in Table El has been omitted. 

Table 2: Values of x 2 /d.o.f. obtained in fits as specified in Table^but with the additional constraint 
D u = D d . 
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We have finally performed a series of fits where we imposed C u = Cd with the constraint C q < 
1.22 GeV -2 , summarized in Table ITBl Provided that we allow a flavor dependence of either f3 q or D q , 
a reasonably good description of the data can be achieved, although the x 2 f° r the neutron data is 
rather high. 

We note that positivity restricts and E% in a rather asymmetric way. The bound (|64[) is 
stronger for d quarks, given the stronger decrease of d v (x) with x. With the profile functions f q (x) 
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Figure 21: Average shift s q {x) of the distance between struck quark and spectators in a polarized 
proton, defined in (fo"4"|) . Solid curves correspond to the default fit (fTT|) . dashed curves to (3 U = 3.5,4.5 
or Pd = 5, 6, 7, and dotted curves to (3 U = 5, 5.5 or (3d = 8, 9. In all cases a = 0.55. 



from our default fit of H$ one needs at least (3d > 5, whereas (3 U may be as small as 3.5. For the 
parameters of the profile function, the constraint (|H5|) is in turn more restrictive on u quarks. We 
have seen that even with these constraints the flavor structure of the forward limit and of the profile 
function in cannot be uniquely determined by present data on the Pauli form factors. The neutron 
data can exclude a complete flavor independence of the distributions (up to their normalization) and 
somewhat disfavors a flavor independent profile function towards larger x. Flavor independence of g„ 
at small to moderate x cannot be excluded from the data, but it requires a rather extreme behavior 
of In view of these uncertainties, precise data on in a wider t range would be of great help. 

In Fig. we plot the average sideways shift s q (x) of the distance between struck quark and 
spectators in a transversely polarized proton, obtained with different values of f3 u and (3d for which 
satisfactory fits to the Pauli form factors can be obtained. Despite the considerable spread of possi- 
bilities, comparison with Fig. [5] shows that | is clearly below d q {x) for all x where we can trust 
our results. 

6 Results for valence GPDs and their moments 

In this section we present the results for the valence GPDs at zero skewness which follow from our 
default fits to the proton and neutron form factors. 

6.1 GPDs as a function of x and t 

In Fig.[22]we have plotted H%(x, t) and H^(x, t) as functions of x for some fixed values of momentum 
transfer t. We recall that the distributions in our default fits refer to the scale /i = 2 GeV. Notice 
the qualitative change when going from small to large values of \t\. For momentum transfer well 
below 1 GeV 2 the GPDs resemble the forward parton distributions and in particular show a divergent 
behavior at small x. For increasing momentum transfer, the GPDs become narrower and develop a 
pronounced maximum, which shifts towards higher values of x in accordance with the behavior of 
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(x)t in Fig. El The height of this maximum decreases with \t\, and this effect is more dramatic for d 
quarks than for u quarks. This can be traced back to our fit result for the profile functions, where 
fd(x,t) > f u (x,t) at large x. The dashed lines in Fig. [221 indicate the values of x where according to 
Fig. [21 the sensitivity of the form factor fits is not sufficient to constrain the GPDs, and therefore the 
result should be viewed as an extrapolation which follows from our particular ansatz for the functional 
form of the x and t dependence. 

The same exercise is repeated for H^(x,t) and H!^(x,t) in Fig. 1231 Since we take the same 
profile functions for H% and H$, the qualitative behavior with increasing \t\ is similar for both cases. 
Remember that for H v our imprecise knowledge of the polarized parton densities represents a major 
source of uncertainty. 

In Fig. [231 we show the corresponding plots for the helicity-flip distributions E%(x,t) and E^(x,t) 
obtained in our default fit to the Pauli form factors and F^ ■ Again we observe a remarkably 
different behavior for u and d quarks. At t = the distributions e%(x) and — e^(x) are not too 
different, but for intermediate \t\ the distribution for d quark develops a maximum which is more 
pronounced and located at significantly smaller x than in the u quark distribution. At larger values 
of \t\ we observe again a faster decrease of the d quark distribution. This is because our default fits 
have equal values for the parameters A q and C q , which respectively govern the large- x behavior of 
the profile functions in H$ and E%. 

6.2 Moments of GPDs 

Important quantities obtained from GPDs are their moments in x. For zero skewness we define 
moments of valence GPDs as 





f 1 dx x n ~ 
Jo 


- l H«{x,t) 


K fl {t) = 


f dxx n ~ 
Jo 






f dxx n ~ 
Jo 


~ l E q {x,t) 



In terms of the notation for moments used in |631 119) we have 

h 9 nfi (t) = A q n>0 (t) , e%{t) = B%{t) for odd n, 

<oW=<oW for even n. (73) 

For even n the moments h? n and are not form factors of local operators, but can rather be seen 
as the valence contributions to the corresponding moments A q n and , which involve both valence 
and sea quarks. Likewise, h q n is the valence contribution to A q n when n is odd. Note that for the 
lowest moments we just have h q = F q and e q = ■ 

In Fig. [23 we have plotted the first three moments of and H^. The scaled u quark moments 
t 2 hno(t) show a similar behavior as the Dirac form factor of the proton, with a smooth increase up 
to values of about —t = 5GeV 2 and a rather flat plateau between lOGeV 2 < — t < 30GeV 2 . On the 
other hand, the lowest d quark moments t 2 h°[ (t) have a pronounced maximum around —t = 3GeV 2 
and die out much faster for higher values of —t. This can be understood in terms of the soft Feynman 
mechanism, or more precisely of the Drell-Yan relation, which relates a large-x behavior like (1—x)" of 
a forward parton distribution with a large-t behavior like \t\~^ 1+ ^i)/ 2 Q f the associated form factor, see 
Sect. 13.51 Since d v (x) behaves approximately like (1 — x) 5 at large x, we expect that t 2 /i^ (0 ~ 
at intermediate values of \t\, which is roughly what we observe in Fig. 1251 
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Figure 22: Result for the valence GPDs H$(x,t) at fj, = 2GeV obtained in default fit to Ff (t) and 
Dashed lines indicate the regions where x < x m ™ (t) or x > x max (t), see (|2*i)) and Fig. [21 
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Figure 24: The same as Fig. |22]fbr the valence GPDs E$. Dashed lines indicate the regions where 
x < x m i n (t) or x > x max (t) according to Fig. IT8l 
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Figure 25: The first three moments of the valence GPDs iJ" (left) and H$ (right), scaled with t 2 . 
The error bands denote the parametric uncertainty resulting from the fit to the Dirac form factors 
Ft and FV. 




Figure 26: Systematic uncertainties for n = 2 moments of valence GPDs H%. The different line styles 
represent alternative fits to the Dirac form factors, as specified in the caption of Fig. 1111 
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-t [GeV 2 ] - t [GeV 2 ] 

Figure 27: The first three moments of the valence GPDs (left) and H$ (right), scaled with t 2 . 
The error bands denote the parametric uncertainty of our default fit. 




Figure 28: The first three moments of the valence GPDs E% (left) and (right), scaled with t 2 . 
The error bands denote the parametric uncertainty resulting from our fit to the Pauli form factors 
iff and F 2 n . 
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The parametric errors on o are rather small, reflecting the quality of the experimental data on 
the proton form factor. On the other hand the uncertainties on /i^ become large for \t\ above, say, 
10 GeV 2 , which is due to the lack of data on the neutron form factor in that region. In addition, 
we have to be aware of the systematic effects related to different choices or constraints for the profile 
functions. In Fig. EH1 we study these effects by comparing the results for the moments h% q from 
different fits to the Dirac form factor data, as we did for the distance function d q (x) in Fig. 1111 We 
see that the systematic error for u quark moments is reasonably small, whereas for d quarks the 
uncertainty on individual moments above \t\ = 3 GeV 2 becomes very large. We should mention that 
in ratios of moments for the same quark flavor the systematic effects drop out to some extent. 

For he and e q n , which are shown in Figs. 1271 and l2"51 similar comments concerning the t depen- 
dence apply. We should point out that the moments h n> o are subject to additional uncertainties due 
to the polarized parton densities Aq v (x), which are not shown for simplicity. In the case of the 
rather large systematic uncertainties discussed in S ect ion 15 . 2 1 should be kept in mind. We remark that 
for our default fit to the Pauli form factors, the powers governing the large- x behavior of the forward 
distributions &%{x) are j3 u ~ 4 and (3d ~ 5.6. Via the Drell-Yan relation this translates into a power 
behavior ~ |t| -2 ' 5 and ~ |i| -3 ' 3 in the region where the soft Feynman mechanism applies. As 
we can see from Fig. 1281 this behavior sets in for values of \t\ above 5 GeV 2 , which is at the border of 
the region presently covered by experiment. Measurements of F% and F% up to momentum transfers 
of, say, 10 GeV 2 could rather directly probe the region of the Feynman mechanism and thus provide 
valuable constraints on the exponents U and fid- The t dependence of the presently available data 
probes the distributions in the transition region between small and large x. 

The lowest moments of GPDs have been calculated in lattice QCD. In particular, the moments 
A?, an d -Sfo f° r u an d d quarks have been obtained in [5] for —t up to 3 GeV 2 . The result has 
been parameterized in terms of a simple dipole fit with a common dipole mass M = 1.11 ± 0.20 GeV 
for all n = 2 moments. One should keep in mind that the value of M has been obtained by linear 
extrapolation to the chiral limit of results obtained with rather high pion (or equivalently quark) 
masses. Furthermore, we have evaluated the valence quark contributions to A^q and -Bf o> whereas 
the lattice calculation includes sea quark effects to the extent that they are included in connected 
diagrams and in the quenched approximation of QCD. In Fig. EH1 we compare our result for the 
normalized moments h\ (*)/^l o(^) an< ^ e 2 o(*)/ e 2 o(^) with the dipole fit to the lattice results. Within 
their uncertainties the results are in remarkable agreement. As already pointed out, our analysis leads 
to quite different behavior of u and d quarks. In particular, q falls of significantly faster than K^q- 
It will be interesting to see whether this effect is also observed in the improved lattice calculations 
currently under way. 



6.3 Valence contribution to Ji's sum rule 

Our determination of the GPDs E$(x, t) from the fit to the Pauli form factors of proton and neutron 
has enabled us to estimate the forward distributions e%{x), which are not accessible in inclusive 
processes but play an essential role in understanding how the total spin of the nucleon is made up 
from quarks and gluons. The forward distributions e%(x) enter the angular momentum sum rule 
in the form 

2(Ll)= ( dx\xe%{x) + xq v {x) - Ag„(x)l , (74) 
Jo 1 J 

where we have only given the contribution from valence quarks. The second moments of unpolarized 
parton densities q v (x) are well-determined from DIS and other inclusive processes. The first moments 
of polarized parton densities can be obtained from the axial- vector couplings of nucleons and hyperons 
under the assumption that the couplings satisfy flavor SU(3) symmetry and that the polarized quark 
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Figure 29: Comparison between lattice QCD results [H] and our default fit for the normalized moments 
h\ (t)//i-2 o(0) an d z\ o(^)/ e l o(0)- The dashed line with the band represents a fit to the lattice data 
with a common dipole mass M = 1.11 ± 0.20 GeV for all n = 2 moments. 

sea is also SU(3) symmetric. For a recent analysis we refer to jSl], whose results were also taken as 
an input for the polarized parton distributions of |21j used in this work. The contribution to (|74j) 
from the second moments of e1{x) can be computed from our default fit. For this particular quantity 
the systematic uncertainties are dominated by our insufficient knowledge on the power (3 q used in the 
ansatz (|67|) . In comparison to this parameter, the power a in (|67|) is rather well determined. 

In Fig. |301 we study the variation of the second moments e\ at t = within the wide range of (3 q 
values for which we found consistent fits to the Pauli form factors. We observe that the contributions 
of u and d quarks are almost equal in magnitude (varying between 0.10 and 0.15) and opposite in 
sign. The corresponding relative uncertainties on the valence contributions to the orbital angular 
momentum are rather small. This is shown in Fig. 1311 where we have used the unpolarized CTEQ6M 
distributions at /i = 2 GeV and the results of to evaluate the last two terms in ()74|) . In the sum 
(L v ) = (L") + (L%) the contributions from u and d quarks cancel to a large extent, so that the relative 
error on this quantity is bigger again. With the central values of (3 U and (3d from our default fit we 
obtain 2{L V ) ps -0.17 at /i = 2 GeV. 

As we already mentioned, sea quarks do not enter the Pauli form factor and are thus not accessible 
in our fit for the GPDs. For completeness we note that the sea quark contributions to 2(L U ) and 2(L d ) 
at fj, = 2 GeV are 2 / dx {xu — Au) ps —0.09 and 2 / dx (xd — Ad) w —0.08 if we take the distributions 
from CTEQ6M [20J and the NLO distributions in scenario 1 of [22 • Since we have no information on 
the sea quark contributions corresponding to the first term in (|74|). we refrain from comparing our 
results with evaluations of (L q ) in lattice QCD [51 l9l I65j. 
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At Ai 
Figure 30: Variation with (3 q of the second moment of e%{x) at \x = 2GeV. 





Figure 31: As in Fig. l3*Ulbut for the valence contribution to 2(L q ) at \i = 2GeV. 
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6.4 Proton tomography 



In Figs. E21 and 0S1 we illustrate our default results for GPDs as tomography plots in impact parameter 
space for fixed longitudinal momentum fraction x. We show the unpolarized density q v (x,b) as well 
as its analog q„(x,b) for a proton polarized in the x direction, see (JHJ and (|50j) . Notice that in 
general q^{x,b) has no rotational symmetry in the impact parameter plane, as is readily seen from 
(|50j) . One can nicely see that the displacement of the center of the distribution along the b y axis is 
different for u and d quarks. According to our discussion in Section El the average displacement is 
described by the quantity s q (x) defined in ()54j) . Indeed we read off from Fig. [^that for small values 
of x the shift |sd(£)| in our default fit is significantly larger than s u (x). This explains the pronounced 
difference between the polarized and unpolarized densities (x, b) and d v (x, b) at small x compared 
to the corresponding u quark densities. The different signs of s u (x) and Sd(x), corresponding to the 
different signs of the anomalous magnetic moments k u and Kj, imply that for d quarks the center of 
the polarized density is shifted toward negative by, whereas for u quarks it is found at positive b y . 

7 Handbag approach to wide-angle Compton scattering 

Our analysis in Sect. 13 .51 has shown that at momentum transfer —t around lOGeV 2 the proton Dirac 
form factor can be described in terms of the soft Feynman mechanism. Dominance of this mechanism 
is the basis of the so-called handbag approach to wide-angle Compton scattering [IJ and related 
processes |661 167j . In this approach the relevant process amplitudes factorize into a hard scattering 
on a single quark in the target and moments of GPDs at skewness £ = 0. The handbag approach 
is restricted to kinematical situations where all Mandelstam variables s, t, u are sufficiently large 
compared to a hadronic scale A 2 , but not yet in the asymptotic regime where the standard hard- 
scattering approach [IT] is applicable. Before turning to a quantitative analysis of the Compton 
form factors on the basis of our results from Section 03 we wish to clarify some important theoretical 
aspects, concerning in particular the inclusion and interpretation of radiative QCD corrections in the 
handbag formalism. A formal proof for handbag factorization would require a systematic treatment 
of these corrections, an issue which has not been addressed in too much detail yet. As a first step, 
one-loop radiative corrections to the partonic scattering amplitude in Compton scattering have been 
evaluated in [H^lfi!?]. 

7.1 Compton scattering and Compton form factors 

Let us recall some important features of the handbag approximation to Compton scattering, jp — ► jp, 
at large angles. (The discussion in this and the following subsection can be generalized to the case 
where the incoming photon is off-shell.) The approximation assumes that the soft configurations 
discussed in Sect. 13.51 are dominant in the scattering process. The process amplitude can then be 
written in terms of hard-scattering amplitudes for Compton scattering on a free quark, jq — > jq, 
multiplied with Compton form factors describing the emission and reabsorption of the struck quark 
by the proton target. The form factors are given in terms of GPDs as 



with the momentum fraction x defined in a frame where the plus-momentum of the proton is equal 
before and after the scattering, see Sect. 13.51 and The derivation of the factorization formula 




(75) 
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Figure 32: Tomography plots of u v {x,b) (left) and u^{x,b) (right) in the transverse b x -b y plane. 
Note that the scale of intensity for longitudinal momentum fraction x = 0.6 differs from the one for 
x = 0.3 and x = 0.05. 
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Figure 33: Tomography plots of d v (x,b) (left) and d^(x,b) (right) in the transverse b x -b y plane. 
Note that the scale of intensity is different for all three longitudinal momentum fractions x. 
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requires approximations appropriate for the soft kinematics specified in (|35|). where terms suppressed 
by A/-y/Ji[ are neglected. 1 In particular, the amplitudes for the partonic subprocess ^q — ► 73 are 
evaluated with the four-momentum k of the struck quark approximated by the full momentum p of 
the proton (see Fig. 112b for the kinematics). The Mandelstam variables s, t, u of the subprocess 
7<7 1Q are then equal to the Mandelstam variables s, t, u of the overall process 775 — ► jp. The 
longitudinal momentum fraction x of the struck quark is thus set to 1 in the hard-scattering amplitude, 
leaving the integration over x to be performed at the level of the proton matrix elements as given 
in (|75|) . The results of our phenomenological fit for H$(x,t) show that at —t of several GeV 2 the 
proton matrix elements are not dominated by very large values of x. In Fig. [22] we see for instance 
that at -t = 10 GeV 2 the maximum of H^(x, t) is around x ~ 0.62, with half the maximum attained 
at x ~ 0.47 and x ~ 0.77. This is well within the parametric estimate (1 — x) ~ A/y^jiy if we take for 
A a value around 1 GeV, but one may worry that the resulting approximation of the hard-scattering 
amplitude is very bad. This need not necessarily be the case. To see this, let us consider more 
sophisticated relations between s, t, u and s, t, u. The ansatz s = xs and u = xu for instance can be 
argued to give a better first approximation than the naive relation with x = 1. Since the leading-order 
scattering kernels for 75 — > jq depend only on u/s, the value of x cancels in this approximation. 2 
A consistent treatment of the kinematics in the handbag graphs beyond the x = 1 approximation 
must however be complemented by the inclusion of graphs that are not described by the proton 
matrix elements corresponding to the twist- two GPDs in (|75|) . Such graphs involve for instance 
operators with two quark fields and an additional gluon, or four-quark operators associated with the 
so-called cat's ears diagrams. The inclusion of higher-twist operators is already required to preserve 
electromagnetic gauge invariance of the Compton amplitude, as has been seen in the treatment of 
deeply virtual Compton scattering [711 172j . A corresponding analysis for wide-angle scattering would 
be very important but has remained elusive so far. 

A related theoretical uncertainty concerns the factor 1/x in the integrals (|75j) . In the derivation 
of [3] it is somewhat ambiguous whether to associate this factor with the soft matrix elements or with 
the hard-scattering amplitude, where it would be approximated by 1. The difference between 1/x 
and 1 is among the effects which are presently beyond theory control. Numerically it is quite large in 
the t range of Compton scattering experiments, given that the expansion parameter A/y / |t| relevant 
for the soft Feynman mechanism only scales with the square root of the large momentum transfer. 
The expressions (j75[) should only be used for values of t that are large enough for the underlying 
approximations — at small t the corresponding integrals even diverge. As we observed at the end of 
Sect. 13.61 the valence GPDs from our ansatz with an exponential t dependence vanish like a power 
of x for x — > provided that \t\ > 0.7 GeV 2 , so that their 1/x integrals then converge. However, sea 
quarks also contribute to the Compton form factors, which describe exchange of positive C parity in 
the t-channel. Since in the distributions at t = the small- x rise is faster for sea quarks than for 
valence quarks, one expects that larger |t| is necessary to make ()75(l well defined. We shall return to 
the issue of sea quarks in Sect. 17.31 

The 1/x factors in (|75|) . as well as the sgn(x) in the case of Ra, give rise to a ^ dependence 
of the Compton form factors. Arising from DGLAP evolution of the GPDs under the integral, this 
scale dependence is an effect of order a s . We shall see in the next subsection that at present NLO 
corrections to wide-angle Compton scattering are only understood in the approximation where x is 
set to 1 in the hard-scattering amplitude and where the form factors (j75fl are replaced with their scale 

1 We recall that the power counting given in [S] refers to the ultrasoft region in 1351 and should be amended accordingly. 
This affects the parametric estimate of corrections to the approximations made, but not their final result. 

2 We remark that such a cancellation does not take place with the expressions for s,t,u recently used in I7U| . These 
expressions cannot be used for wide-angle Compton scattering, since they imply a zero of u at a particular value of x, 
see Sect. 4.1 of [5]. 
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independent counterparts 



RvV = E e * f dxHi(x,t,v 2 ) , 4 0) (*) = E e * f dxEi(x,t, 

q *' — 1 q J -1 

4° } W = E e * f_teH*{x,t,f), (76) 

which apart from charge factors are just the usual Dirac, Pauli, and axial form factors. In addition 
to neglecting the difference between 1/x and 1, we have counted antiquarks with the "wrong" sign 
in the integral over H q , which should again be a valid approximation when large enough momentum 
fractions dominate in the integral. 

7.2 NLO corrections to wide-angle Compton scattering 

Next-to-leading order corrections to the partonic subprocess in Compton scattering have been calcu- 
lated in jSH] • This means to consider virtual corrections to jq — > jq where the virtualities of the loop 
momenta are hard, i.e., where they scale with the large momentum transfer t. The unrenormalized 
result for the hard-scattering amplitudes has the schematic form 3 

n LO (s,t) (i + ^ Qr(^ 2 )) + ^ w NLO M) , (77) 

where the NLO contribution has been decomposed into a divergent part Cir(/U, t) multiplied with the 
LO amplitude and a finite piece W NLO (s,t). Using dimensional regularization with D = 4 + e, one 
has 

ClR,(i, fj 2 ) = — j H — In — H In 2 — — 3 In — + constant terms , (78) 

which contains a collinear and a soft infrared divergence, and the well-known Sudakov double log- 
arithms. It has been pointed out in [SH] that the soft and collinear divergences and the associated 
double logarithms are universal: not only do they appear in the corrections to 'yq — ► 'yq but also in 
the corrections to -y*q — ► q, which are relevant for the elastic proton form factors. In a "physical" 
factorization scheme one considers the radiative corrections to the quark-photon vertex as part of the 
form factor itself. To renormalize the Compton form factors and the corresponding hard-scattering 
amplitudes one then subtracts the full expression Cir, where the constant terms are fixed by the con- 
dition that in this scheme the elastic proton form factor is unchanged. 4 Writing the overall Compton 
amplitudes as Mi, where i corresponds to the subscripts V, A, T of the Compton form factors, we 
then have 

/L(V +\ i a sCF -tjNLO/^ + \\ r> (a ,,2\ 



Mi(s,t) = (n^( s ,t) + ^n^ u (s,t)jR i (t,^) + ... 

nf%, t) + ^ nf^(s, *)) itf (t) + nl°( s , t) [r, (t, - R f (t)} +.... (79) 

Note that the NLO calculations for 'yq — > -yq and 7*g — > q with on-shell quarks enabled us to find 
a renormalization scheme at the level of the electromagnetic and weak nucleon form factors and of 
their flavor combinations (|76|) relevant for Compton scattering. They do not permit an analogous 

3 Our discussion refers to the case of photon helicity conserving amplitudes. Photon helicity-flip amplitudes are zero 
at LO and finite at NLO. 

4 The relation between our notation in and the one in EH1 is Hj^ = - (n^t° + it - 8 



51 



renormalization of the Compton form factors (|75j) or of the GPDs themselves. This is manifest in 
the first line of (|79|). where we have an uncanceled yi dependence. In the second line of (|T§|) we 
have therefore indicated that, for internal consistency of the presently available results, the NLO 
calculation should be restricted to the part of the Compton form factors which involves the same 
local operators as F±, F2 and Fa- The parametric uncertainties of the second line in (|79|) are of order 
a 2 or of order a s A/y/\t\. Up to this accuracy the result is 11 independent as it must be. Controlling 
effects of order a s A/>/|i[ would require a more complete analysis of power-suppressed terms already 
at Born level. 

Following an alternative point of view, one could use a minimal subtraction scheme in (|77|) to 
render the resummation of Sudakov logarithms explicit in hard coefficient functions. The Compton 
amplitudes would then be written as 



M i {s,t) = C i (s,t,Li A )Rf (t,^) + ..., (80) 
where Ci(s,t, li 2 ) obeys the renormalization group equation 



V 2 ^2 C ^ s i *- = 7(t, ^) C t (s, t, li 2 ) (81) 



d/i 

whose anomalous dimension follows from (|78|) 




7 (*,M 2 )=-^ In^r+Tb \+0(a 2 ). (82) 



The value of the constant 70 depends on the renormalization scheme. Renormalization group evolution 
from 11 2 ~ |t| down to ix\ ~ Ay^jtj corresponds to the summation of Sudakov logarithms between those 
two scales, which are then explicitly contained in d(s,t, liq). In this sense, the new Compton form 
factors Ri S (t, liq) only contain dynamics below the factorization scale liq, which is the virtuality of 
the struck quark in soft kinematics. They can be considered as being defined in soft-collinear effective 
theory [731 174j , which is an effective field theory where the hard degrees of freedom are integrated 
out, leaving degrees of freedom with virtualities smaller than \t\. The scale dependence of i?| ff (t, // 2 ) 
compensates that of Ci(s,t, li 2 ). 

Similarly, in the effective theory the Dirac form factor can be written as 

F 1 (t) = C 1 (t,Li 2 )Ff(t, L i 2 ) + ... (83) 

with a hard coefficient function that obeys the same evolution equation as in the case of Compton 
form factors, 

C ^ ^) = 7& ^) C ^ M 2 )- ( 84 ) 

In other words, the function C\{t, Li 2 ) resums the same (leading) Sudakov logarithms as Cj(i, li 2 ). The 
physical scheme where the coefficient function in (|83[) is equal to 1 is related to a general scheme by 
a finite renormalization. Applying the same renormalization to the Compton form factors one would 
write 

Rf ) \t) = C 1 (t,» 2 )Rf(t,n 2 ) + ..., (85) 

which together with (|8U|) gives 



Cj(s,t,Ll 2 ) D (o), 
Ci{t,ii 2 



M l (s,t) = ^;;\,J Rr(t) + -..- (86) 
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Sudakov effects cancel in the ratio of hard-coefficient functions, 5 which can therefore be calculated in 
fixed-order perturbation theory as 

C^y?) ~ Hi ( S '*) + ^T^ + (87) 

where we recover (|79[) to leading accuracy in the expansion parameter A / y/\t\ of the effective theory. 

Returning to the level of GPDs, we have so far dealt with distributions H$(x, t, /x 2 ) defined by 
standard collinear renormalization in QCD. Note that even for /i 2 <C \t\ such distributions can be 
defined, since they can be obtained from distributions with /j 2 > \t\ by DGLAP evolution, which 
may be understood as a finite renormalization. The result does however not correspond to matrix 
elements with internal virtualities of at most A-y/jtj", which would be appropriate for the soft Feynman 
mechanism. One may speculate that such a "soft GPD" would be obtained in the effective theory as 

H**{x,t,n 2 ) = _l_fl?(x,t,^) (88) 
for \t\ 3> A 2 , with associated sum rules 

t Jo t Jo 

Notice that H$ eS obeys a different evolution equation than H$, since its renormalization includes 
the Sudakov logarithms from the endpoint region x — ► 1. It is however not clear whether the x- 
independent renormalization factor in (|88|) corresponds indeed to the structure emerging in the effec- 
tive theory. 

In the absence of a proper definition for "soft GPDs" , we understand the distributions extracted 
in this work as "standard GPDs" obeying DGLAP evolution. This is also convenient in connection 
with our phenomenological ansatz (jllj) . which interpolates between the standard MS parton densities 
at t = and distributions at large t. We point out that even for the highest values of —t ~ 30 GeV 2 
we consider in our applications, Sudakov double logarithms are not very large, with ln 2 (/x 2 /|t|) < 4 
for (j, = 2 GeV, so that one may expect the difference between "soft" and "standard" GPDs to be 
reasonably small. 

In our calculations for wide-angle Compton scattering we will use the form factors (|75|). regarding 
their 1/x factors as a phenomenological estimate of effects beyond a strict A/ y/\t\ expansion. To keep 
the analysis simple we have refrained from using the expression ()79|) and instead take the full form 
factors (|75|) multiplied with the hard-scattering amplitudes at NLO. Since the /j, 2 dependence inherent 
in the 1/x moments is an effect beyond present theoretical control, we will require that the variation 
of the form factors obtained with GPDs fitted for /x = 1, 2 and 4 GeV remains reasonably small. This 
will provide a criterion of how large t is required to obtain stable results within our present approach. 



7.3 Results for form factors and observables 

After our discussion of theoretical uncertainties let us now see how they quantitatively look like with 
the GPDs we have obtained from our fits to elastic form factors. We restrict this discussion to Ry, 
which dominates the unpolarized Compton cross section. In Fig. we compare the full form factor 
Ry and its analog Ry^ without a factor of 1/x under the integral, both evaluated with the result of 

The situation is similar to the case of heavy-to-light form factors, where the issues related to the summation of 
Sudakov logarithms in the effective theory are irrelevant for ratios of soft form factors |75|. 



53 



1.5 




°0 2 4 6 8 10 12 14 
-t [GcV 4 ] 



Figure 34: Study of systematic uncertainties regarding the Compton form factor Ry. Shown are 
t 2 Ry and t 2 Ry^ from (|75|) and (|76|) evaluated from our default fit for H$ at \i = 2 GeV. The limits of 
the error bands correspond to the fits for H% at fx = 1 and 4 GeV in Table I1UI The dotted line shows 
the sea quark contribution to t 2 Ry as specified before and in (|9U|). 



our default fit for H$(x,t). As is to be expected from the shape of H$ in Fig. [22 the difference is 
considerable in the t range relevant for existing or planned measurements. The error bands on the 
curves for Ry and Ry ^ indicate the variation of H$ with fi, with the limits of the bands corresponding 
to the fits with \i = 1 and 4 GeV listed in Table 1101 The scale variation for Ry becomes rather large 
for —t below about 3 GeV 2 , which gives a first indication of when effects of order a s K/^/\t\ can no 
longer be neglected. The tiny scale dependence of R^y confirms that our fits to Ff and F™ with input 
parton densities at different fx indeed give a scale invariant result to a very good approximation. 

As we already mentioned, sea quark contributions do not cancel in the Compton form factors 
defined in (|75|). Contrary to the isovector axial form factor Fjy(t), the quark flavor combination 
in Compton scattering is also sensitive to the flavor singlet distribution, which is more difficult to 
model due to its mixing with gluons as explained at the end of Sect. 13.11 For a naive estimate we 
have taken the same impact parameter profile for sea quarks and valence quarks. We thus define 
HQ(x,t) = -H q (-x,t) and set H^(x,t) = q(x)exp[tf q (x)] with the CTEQ6M j2U] antiquark densities 
q(x) at /i = 2 GeV and with f q (x) determined in our default fit for H%. The resulting sea quark 
contribution 

i?r(i,/?) = 2]Te 2 F^W&t,^) (90) 

q JO X 

to Ry is shown as the dotted curve in Fig. 1341 where we have restricted the sum over quark flavors u 
and d. When the sea quark contribution becomes important, we are clearly outside the region where 
the physical picture and the approximations of the soft Feynman mechanism are applicable. We see 
that with our estimate one need not worry about sea quarks for —t > 4 GeV 2 . Below —t = 2 GeV 2 
their contribution quickly grows out of control, and we will therefore not show Compton form factors 
below this value. 

We have seen in Fig. El that our default fit for H$ exhibits the scaling properties characteristic of 
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the soft Feynman mechanism in Ff(t) for —t above 10 GeV 2 . Below this value one enters a transition 
region, where the scaling becomes more and more approximate. Figure |3U indicates that current 
theoretical control over the Feynman mechanism in wide-angle Compton scattering does not extend 
below -t = 2 GeV 2 . 

In Fig. OH] we show the Compton form factors obtained from our default fits for E% and from 
the ansatz (|46|) for H$ with f q = f q . We plot both individual flavor contributions and the flavor 
combinations R{ = ^Rf + ^Rf for Compton scattering on the proton. In the case of R v and R q A 
we see the clear dominance of u quarks over d quarks even without the charge factors, as we already 
did for the corresponding moments in Sect. 16.21 The full form factors Ry and Ra turn out to be 
somewhat smaller than those modeled in (Hj. 

For Rj, the relative weight of flavors is different than for Ry and R q A . Even more than for the 
moments ef in Fig. 1281 d quarks dominate for smaller and u quarks for larger values of t. Given the 
opposite signs of R T and R T and the squared charge factors, the resulting Rt is significantly smaller 
than Ry and Ra at lower values of t. The ratio Rt/Rv behaves quite differently from the analogous 
ratio F^/Ff of electromagnetic form factors. With 

Kc{t) = ^ rW) ' KemW = l^ifw' (91) 

the measurement of [JHl 

gives an approximately constant K em ra 0.37 for 2 GeV 2 < -t < 5.5 GeV 2 , 
whereas we find that with our default fits k,q rises from 0.1 to 0.25 in the same interval and continues 
to grow to 0.4 at —t = 10 GeV 2 . Our result for Rt is however subject to rather large uncertainties, 
given the considerable freedom we encountered in extracting E$ from the Pauli form factors alone. 
The fits shown in Table WH which all give a good description of F% and i 7 ^, produce form factors Rt 
which at —t = 5 GeV 2 are up to a factor of 1.5 larger than Rt shown in Fig. I35( with even larger 
discrepancies at lower t. Within this set of fits, the curve in the figure acts more like a lower bound. 

With the form factors in Fig. OH] we obtain the unpolarized cross section da/dt shown in Fig. 1361 
for a squared cm. energy of s = 11 GeV 2 . We have included the NLO corrections in the physical 
scheme described in Sect. 17.21 subtracting a constant term as ^ F — 8^ from the quark scattering 
kernels of [HE] as mentioned in footnote B] For the form factor Ry describing gluon exchange at NLO 
we have taken the model result of [HE] divided by a factor of 1.3, which is the typical discrepancy at 
— t ~ 5 GeV 2 between Ry obtained in the present study and in |68| . The resulting gluonic contribution 
to da/dt is only a few percent. 

The inner band for the curve of da/dt reflects the parametric errors on the form factors shown in 
Fig. 1351 The corresponding uncertainty is very small. This is because Ry is strongly dominated by 
H%, which is the distribution best constrained by the abundant data for Ff, whereas the individual 
contributions of both Ra and Rt to da/dt are between 1% and 10% for the kinematics shown in 
Fig. 1361 The central curve in the figure corresponds to scenario 2 in |77j . where the Mandelstam 
variables of jq — > jq and of jp — ► 7p are related by s = s — m 2 , t = t and u = u — m 2 . The limits of 
the outer band correspond to scenarios 1 and 3 in the same study and reflect the uncertainty in the 
evaluation of da/dt that is due to the finite proton mass. 

An observable with greater sensitivity to Ra is the correlation parameter All between the he- 
licities of the incoming photon and the incoming proton [781 168j , which to a good approximation is 
given by the ratio Ra/Rv times a known kinematical factor. Within our approach All is equal to 
the correlation parameter Kll between the helicities of the incoming photon and the outgoing pro- 
ton, which is a consequence of neglecting the quark mass in the hard-scattering subprocess (see also 
|79j). Our result for s = 11 GeV 2 is shown in Fig. 1361 with the error band reflecting the parametric 
uncertainties on Ra, Ry and Rt- 
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Figure 35: Scaled Compton form factors for individual quark flavors and summed with the appropriate 
squared charge factors (bottom right). The error bands shown for Ry and for Rt correspond to the 1 a 
uncertainties of our default fits for H% and E$. The band for Ra corresponds to the 1 a uncertainties 
on our profile function added in quadrature to the errors on the polarized parton densities in |21j . 
which are by far dominant for this quantity. 
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Figure 36: The unpolarized cross section (left) and the helicity correlation parameter All (right) 
for wide-angle Compton scattering at s = 11 GeV 2 as a function of the scattering angle 9 in the cm. 
Both observables are evaluated at NLO with the Compton form factors shown in Fig. 1351 The error 
bands are explained in the text. 

The correlation parameter Kls between the helicity of the incoming photon and the transverse 
polarization of the outgoing proton is sensitive to the tensor form factor Rrp- Assuming kq = n cm ~ 
0.37 the study found values well in the range of 10%. With our default fit we find Kls at 
s = 11 GeV 2 as small as 1% for backward angles 6 and yet smaller in the forward hemisphere. This 
can be traced back to the approximate cancellation of two terms in the expression for Kls, one 
proportional to Rt and the other to Ry, see eq. (397) in ^2]- For this particular kinematics, the 
uncertainties on Rt from our fits hence greatly amplify in the observable Kls, which unlike All 
turns out to be strongly energy dependent. 

8 Summary and outlook 

We have developed a parameterization of valence quark GPDs at zero skewness which takes into 
account dynamically motivated correlations between their x and t dependence. At small x their form 
is consistent with dominance of the leading meson trajectories known from Regge phenomenology. 
For large t they incorporate the soft Feynman mechanism, where the struck quark has a virtuality 
of order Ay^ji[. Both analytically and numerically our distributions satisfy the Drell-Yan relation 
between the large- x power behavior of parton distributions at t = and the large-t power behavior 
of the associated elastic form factors. This holds for both H% and E% and is to be understood in the 
sense of effective powers describing the x or t dependence in a finite range of these variables. Key 
features of our functional ansatz for the GPDs are approximately invariant under DGLAP evolution, 
so that with an appropriate change of the free parameters this ansatz can be used for a reasonable 
range of scales (i. 

Good fits of the GPDs can be obtained to describe both the Dirac and Pauli form factors of 
proton and neutron, with an accuracy better than 20% for almost all data points and better than 5% 
in wide ranges of t. From Ff, where most data is available, we get a good determination of H%(x,t) 
up to x<0.8. It is stable under a change of parameterization, provided we assume an exponential 
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t dependence as in Under this assumption we have thus well determined the average squared 



the distribution H^(x,t) is less well determined, but we find clear indications that (b 2 )^. > (b 2 )" for 
larger values of x. This implies that with increasing \t\ the suppression of d quarks compared with u 
quarks becomes even stronger than in the forward parton densities at large x. 

With the fit results for H% we also achieve a rather good description of the isovector axial form 
factor Fa under the simplifying assumptions that sea quarks can be neglected and that the impact 
parameter distribution of valence quarks is independent of the quark helicity. With our ansatz for 
H% we find that the largest uncertainty on the resulting form factor Fa is due to the current errors 
on the polarized quark densities Ad and Au. 

Using our ansatz for the proton helicity flip distributions we get good fits to the Pauli form factors 
F% and F£. We find however considerable ambiguities in determining E%(x,t) since its forward limit 
is not known. The form factor data disfavor an identical shape in x and t of E%(x,t) and E^(x,t), 
but they have only little preference concerning the details of this flavor dependence. In particular, 
we cannot determine whether the forward limit E„(x,0) has a harder or a softer x dependence 
than E$(x,0). Positivity constraints prove to be extremely valuable in constraining E% and E$, as 
they severely restrict their size for x>0.5 as well as their t dependence. We find that, despite the 
uncertainties on these distributions, their moments / dx xE%(x, t) are reasonably well determined 
within our ansatz. With Ji's angular momentum sum rule, the resulting valence quark contribution 
to the orbital angular momenta (L u ) and (L d ) have even smaller relative errors for the allowed range 
of our parameters. 

The x moments of our extracted GPDs are in reasonable agreement with results from lattice QCD, 
given the uncertainties of both our phenomenological determination and of the lattice calculations. 
With our results we have also calculated the distribution in the transverse plane of valence quarks with 
momentum fraction x, both in an unpolarized and in a transversely polarized proton. We have thus 
shown that "proton tomography" as proposed in |231 124j is feasible on the basis of experimental data. 
Finally, we have used our GPDs to evaluate the form factors needed in the soft handbag approach to 
wide-angle Compton scattering, as well as the associated experimental observables. From a number 
of indicators we conclude that consistency of this approach requires momentum transfers —t above 
3 GeV 2 or so. 

Where can further progress be made? High-quality data on the neutron form factors in a wide t 
range would be highly valuable for pinning down the differences in the spatial distribution of u and d 
quarks. Our fit results indicate in fact drastic differences in the behavior of u and d contributions to 
the form factors, which should be tested experimentally. The fit to the Pauli form factors suggests that 
the t behavior of F^ will continue to change beyond 5.5 GeV 2 , and it would be good to see whether this 
is indeed the case. Since Dirac and Pauli form factors are most natural to connect with the physics of 
parton distributions, we encourage experimenters to provide data on the Sachs form factors Gm and 
Ge at equal values of t and with correlated errors wherever possible. Only then can one transform to 
F\ and F% without introducing further theoretical bias or having to guess errors. With the available 
data on the axial form factor Fa we could only take a glimpse at the role of quark helicity. More 
data on this fundamental quantity (in the form of data points rather than dipole parameterizations 
whenever possible) should lead to important progress in the study of H q . An improved determination 
of the polarized forward densities Aq is under way at several experimental facilities and will also be of 
benefit to this study. For a more precise study of wide-angle Compton scattering, progress is required 
in the theory of power-suppressed contributions. Together with ongoing and planned measurements 
this should lead to an improved understanding of the dynamics and, if the soft handbag approach is 
further validated, eventually to pinning down the relevant form factors. 
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The fits to the electromagnetic form factors in this work have restricted us to zero skewness and 
to the valence quark sector. The ansatz we have developed for GPDs at £ = can be used as an input 
for modeling distributions at nonzero £, using the concept of double distributions as shown e.g. in |80j . 
A nontrivial interplay between t and the longitudinal variables can readily be implemented in this 
framework |81| 118) . Interesting observations concerning this formalism for proton distributions have 
recently been made in [S2], building on earlier work in |831 184) . Most of the exclusive processes where 
GPDs at nonzero £ occur are sensitive to sea quarks and also to gluons, in particular the important 
channels of deeply virtual Compton scattering and of neutral vector meson production. It will be 
interesting to investigate the dynamical interplay between the quark singlet and gluons in the impact 
parameter distributions, similarly to what has been seen for the distributions in momentum fraction 
x. In some channels, like the production of pseudoscalar mesons or pion pairs, one can however 
selectively probe the valence quark combinations of GPDs, see Sect. 1 of |18j and Sect. 5 of [EH- This 
opens the possibility to investigate the transverse distribution of valence quarks at small x, where 
elastic form factors have only a limited reach. The versatility of the GPD formalism and vigorous 
experimental activity let us hope that important progress is still to come. 
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A Details on form factor data 

Experimental data are usually given for the magnetic and electric form factors Ge and Gm, or for 
their ratio. Natural quantities for our parameterization of GPDs and their physical interpretation are 
the Dirac and Pauli form factors, so that we must perform the conversion 

G E + tG m „ G M - G E , . 

Fl = 1 + r ' F2 = 1 + r ' (92) 

where r = —t/(4m 2 ). To obtain central values for F\ and F2 is rather straightforward, but we could 
not derive errors on these quantities in a rigorous way, lacking information on the correlation between 
the errors on Ge and Gm- We encourage experimentalists to provide correlated errors on the two 
form factors (or derived quantities) whenever possible, since depending on the physical context one 
or the other set of form factors is more adequate. 

For the proton form factors we have used the results from [S3], where the original data from 
[12 EE1 E3 was reanalyzed using the additional information on the ratio Ge/Gm measured by the 
recoil polarization method in [881 176) . With increasing t the measurements using this method deviate 
from those obtained with a Rosenbluth separation, which is affected by QED radiative corrections 
growing with t (see |70) and references therein). 

To obtain Ff we have used the central values and errors on G P M together with the analytic 
parameterization for G P E /G P M given in j^SJ. In calculating errors on F\ we have only used the errors 
on Gm but not those on the parameterization of Ge/Gm- The uncertainty on Ge/Gm is included 
in the errors on the reanalyzed Gm data given in [HS], and counting it again in the conversion 1)92)1 



59 



Table 3: Details of the data sets for which separate values of x 2 are given in the tables of Add. IBl 



data set 


references 


data points 


range of —tj GeV 2 


Ff low 




49 


0.156 to 7 


Ff high 


m 


13 


2.862 to 31.2 




[HI EHl ED El EU 


8 


0.255 to 4 


F\ total 




70 




Fl low 


m 


18 


0.0389 to 1.75 


Ff high 


nana 


28 


0.32 to 5.54 




[Ml El EH E3 EU 


8 


0.255 to 4 


F 2 total 




54 





would not be consistent. In discarding them we obtain the same relative errors on Ff as those on 
Gm- They are between 0.8% and 2.4% for —t < 10 GeV 2 and then grow up to 9% at the highest t of 
the data. If we had added the errors on the parameterization of Ge/Gm in quadrature, the resulting 
errors on Ff would have been increased by at most a factor of 1.5 for —t < 7.7 GeV 2 . (Beyond this 
value the ratio Ge/Gm is taken as zero in |85) and it is unclear which error one should assign to it.) 

For Ff we have used the analytic parameterization of G P M given in [S3] together with the data 
and errors on G P E /G P M measured with the recoil polarization method in |88l 1761 159*] . For calculating 
the errors on F 2 we have used those on Ge/Gm but not those on the parameterization of Gm (which 
overestimate the experimental errors on Gm in the t range we need). If we were to add in quadrature 
an error on 

G M of 1.5% for -t < 1 GeV 2 and of 2% otherwise, our errors on Ff would increase by 
at most a factor of 1.6. The recoil data only go down to —t = 0.32 GeV 2 , and in order to have some 
data for lower — t we have also included the results of [0U] for Gm and Ge obtained by a Rosenbluth 
separation, using (|9*2*j) and adding the errors on Gm and Ge in quadrature. The values for Ff we 
obtain from these data match rather well with those from the recoil measurements in the t region 
where they overlap (see Fig. 117(1 . indicating that in the kinematics of interest the two methods are 
reasonably consistent. The errors we estimate for Ff are between 1% and 4%, with the exception of 
two data points with 5% and one point with 7%. 

For the neutron form factors we have restricted ourselves to results published after 1990 and 
only considered values of t where data for both G M and G E are available, adding their errors in 
quadrature (as well as statistical, systematic and theoretical errors when they are given separately). 
The measurements in [611162) cover a range 0.5 GeV 2 < — t < 1.5 GeV 2 , and we obtain Ff with errors 
between 6% and 17% and Ff with errors between 0.8% and 2.4%. At -t = 0.255 GeV 2 we have 
combined the measurement of G M from [""""j] and the one of G E from [HU] . The resulting error on Ff 
is 5%, whereas we obtain Ff = —0.010 ± 0.035 with a huge error. This is because the individual 
terms G% = 0.066 ± 0.037 and tG m = -0.077 ± 0.003 in nearly cancel. We nevertheless keep 
this data point, since a x 2 fit to the data will correctly take into account its large uncertainty. At the 
high-t end we have taken data from j§"f"j. We do not use their results at -t = 1.75 GeV 2 and 2.5 GeV 2 
since the central values of (G E ) 2 came out negative. From the remaining points —t = 3.25 GeV 2 and 
4 GeV 2 we obtained Ff with errors around 30% and F% with errors around 20%. 

To have an indication of how well our fits describe proton data in the lower and upper range of 
available t, we give in App.[B] separate values of x 2 f° r data sets labeled "low" and "high" as specified 
in Table "a 
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B Details of fit results 



In the Tables of this section we list the results of the various fits discussed in this paper. We have 
performed x 2 fits to the form factor data using the minimization package Minuit [H2]- The la errors 
on the parameters are defined by the contour in parameter space for which 

X 2 = xLn + ^X 2 , (93) 

where A% 2 corresponds to a confidence level of 68%. One respectively has A% 2 = 2.30, 3.53, 4.72 for 
2, 3, 4 free parameters in the fit. The error bands shown in our figures have been derived from the 
la errors on the fitted parameters by standard error propagation (see e.g. Sect. 31.2 of [S'2\). This 
takes properly into account correlations between errors on different parameters, which are typically 
quite large in our fits. 

Our default fit for H$ is specified in Sect. 13 .21 and listed in the first row of Table |B1 and our default 
fit for E% is explained in Sect. 15.31 and shown in Table 1131 The covariance and error correlation 
matrices of these fits are given in Tables H] and 03 respectively. All our fits for E% use the exponential 
ansatz (|57|l with the profile function (|61|) and the forward limit (|67|) . Except for Table El we have 
fixed a = 0.55 in all fits. No error is given on C u when a fit has chosen its maximum value 1.22 GeV -2 
allowed by the positivity constraint (|65|) . 

In Tables l6l to 1161 we give the fitted parameters and their la errors. For the sake of legibility we 
have omitted the units of A q , B q , C q , D q and a', which are always taken as GeV -2 . The columns 
labeled "low" and "high" give the x 2 P er number of data points for the separate proton data sets 
specified in Table 01 and the columns "p" and "n" give the corresponding numbers for the combined 
proton data and for the neutron data. The column "total" gives the x 2 P er degrees of freedom, which 
is the quantity minimized in the fit. Note that if different fits in the same table have different numbers 
of free parameters, then the normalization of x 2 /d.o.f . differs slightly from one row to another. The 
precise x 2 values should be interpreted with care because of the inadequacies affecting our errors on 
the Dirac and Pauli form factors (see App. [XJ. This applies in particular to the low-i data on Ff, 
whose errors are particularly small. To assess the quality of a fit it is hence useful to consider in 
addition to its x 2 a l so the- "pull" as shown in Figures [7| and 

From Table 01 it is clear that the low-t data on Ff tend to dominate the x 2 minimization in the fit 
to the Dirac form factor because of the large number of data points. This could be circumvented by 
reweighting the x 2 f° r different data points, but we have opted for the more straightforward procedure 
without weights, given that our best fit turns out to describe the large-t data on Ff and the data on 
Fi quite well. In a similar fashion, the proton data tend to dominate over those for the neutron in 
the fits to both the Dirac and the Pauli form factors due to the paucity of neutron data. 

Table 4: The covariance matrix (left) and correlation matrix (right) of the default fit to H$ described 
in Sect. 13.21 and listed in the first row of Table |H1 The entries of the covariance matrix are in units of 
GeV -4 and are normalized such that diagonal elements give the squared 1 a errors on a parameter. 





B q 


A u 


A d 




B q 


A u 


A d 


B q 


6.95 • 10" 4 


-4.84 • 10" 4 


5.69 • 10~ 3 


B q 


1 


-0.743 


0.750 


A u 


-4.84 • 10~ 4 


6.11 • 10~ 4 


-1.38 • 10~ 3 


A u 


-0.743 


1 


-0.194 


A d 


5.69 • 10~ 3 


-1.38 • 10~ 3 


8.22 • 10~ 2 


A d 


0.750 


-0.194 


1 
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Table 5: The covariance matrix (left) and correlation matrix (right) of the default fit to E% described 
in Sect. 15.31 and listed in Table IT3l The entries of the covariance matrix are in units of GeV~ 4 and 
are normalized such that diagonal elements give the squared 1 a errors on a parameter. 







D u 


D d 




0u 


D u 


D d 




5.01 • 1(T 2 


-2.15 • 10~ 2 


-8.24 • 10~ 3 




1 


-0.884 


-0.693 


D u 


-2.15 • 10~ 2 


1.17- 1(T 2 


1.80 • 10" 3 


D u 


-0.884 


1 


0.311 


D d 


-8.24 • 10~ 3 


1.80 • 10~ 3 


2.84 • 10~ 3 


D d 


-0.693 


0.311 


1 



Table 6: Fits to H$(x, t) = q v {x) exp[tf q (x)] using the CTEQ6M parton densities at scale fj, = 2 GeV 
and f q (x) = a'(l — x) n ~ l log(l/x) + (A q — a')(l — x) n from ()23|) and ()27|) . More information on the 
different columns is given in the text. 



11 


fit parameti 
a' 


3rs 


A d 


X 2 /data points 
low high p 


n 


X 2 /d.o.f. 
total 


1 

2 


1.38 ±0.01 
1.08 ±0.02 


0.070 ± 0.006 
0.78 ±0.01 


0.25 ±0.04 
1.05 ±0.05 


1.42 0.60 1.25 
4.48 8.79 5.39 


1.55 
2.22 


1.34 
5.25 



Table 7: Fits to H^(x, t) = q v (x) exp[tf q (x)] using the CTEQ6M parton densities at [i = 2 GeV and 
f q {x) = q'(1 - x) 2 log(l/x) + B q (l - x) 2 + A q x(l - x) from with fixed a' = 0.9 GeV" 2 . 



fit parameti 
B u = B d 


srs 

A u A d 


X 2 /data points 
low high p 


n 


rVd.o.f. 

total 


0.43 ±0.02 
0.32 ±0.01 


0.164 ±0.014 1.35 ±0.21 
0.190 ±0.009 


3.05 1.75 2.77 
4.23 1.64 3.69 


2.15 
9.46 


2.82 
4.47 



Table 8: Fits to H^(x, t) = q v (x) exp[tf q (x)] using the CTEQ6M parton densities at /i = 2 GeV and 
f q {x) = q'(1 - x) 3 log(l/x) + B q (l - x) 3 + A q x(l - x) 2 from © with fixed a' = 0.9 GeV" 2 . 



fit parameters 




X 2 /data points 




X 2 /d.o.f. 


B u B d 


A-u A d 


low 


high 


p 


n 


total 


0.59 ±0.03 


1.22 ±0.02 2.59 ±0.29 


2.27 


0.70 


1.94 


1.12 


1.93 


0.59 ±0.03 0.32 ±0.14 


1.26 ±0.04 3.82 ±0.78 


1.78 


1.32 


1.68 


0.73 


1.67 


0.56 ±0.02 0.75 ±0.08 


1.17 ±0.03 


3.22 


0.76 


2.70 


2.36 


2.78 


0.50 ±0.01 


1.20 ±0.02 


2.98 


0.69 


2.50 


9.24 


3.36 
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Table 9: Fits to H%(x,t) = q v (x) exp[tf q (x)] with f q (x) as in Tableland different sets of CTEQ6M 
parton distributions at jjL = 2 GeV. 





fit parameters 




X 2 /data points 




X 2 /d.o.f. 


set 


B u = B d 


A u 


A d 


low 


high 


P 


n 


total 


17 


0.55 ±0.03 


1.21 ±0.02 


2.52 ±0.30 


1.79 


1.23 


1.68 


1.09 


1.68 


18 


0.64 ±0.03 


1.24 ±0.02 


2.60 ±0.27 


2.96 


0.58 


2.46 


1.30 


2.44 


35 


0.58 ±0.03 


1.25 ±0.03 


2.88 ±0.28 


2.20 


0.72 


1.89 


0.98 


1.86 


36 


0.59 ±0.03 


1.22 ±0.03 


2.85 ±0.29 


1.99 


0.91 


1.77 


0.81 


1.73 



Table 10: Fits to H$(x,t) = q v (x) exp[tf q (x)\ with f q (x) = a'(l -x) 3 log(l/x) + B q (l -xf +A q x(l - 
x) 2 from (|29j) and CTEQ6M parton densities taken at different scales /i. 





fit parameti 


3rs 






X 2 /data points 




x7d.o.f. 




a' 


B u = Bd 


A u 


A d 


low 


high 


P 


n 


total 


lGeV 


1.17 ±0.05 


0.50 ±0.12 


1.81 ±0.07 


2.62 ±0.31 


1.16 


0.70 


1.06 


1.11 


1.13 


2 GeV 


0.97 ±0.04 


0.43 ±0.09 


1.32 ±0.06 


2.51 ±0.35 


1.77 


1.53 


1.72 


0.86 


1.72 


4 GeV 


0.89 ±0.03 


0.36 ±0.08 


1.10 ±0.06 


2.58 ±0.38 


2.23 


2.09 


2.20 


0.89 


2.17 


8 GeV 


0.84 ± 0.03 


0.29 ±0.07 


0.96 ± 0.06 


2.71 ±0.41 


2.56 


2.42 


2.53 


1.10 


2.51 



Table 11: Fits to fl"«(x,t) = q v (x) (1 - tf q (x)/p)- p with f q (x) as in Table EH and the CTEQ6M 
parton densities at \x = 2 GeV. 





fit parameters 








X 2 /data points 




rVd.o.f. 


p 


a' 


B u = B d 


A u 


A d 


low 


high 


P 


n 


total 




0.97 ±0.04 


0.43 ± 0.09 


1.32 ±0.06 


2.51 ±0.35 


1.77 


1.53 


1.72 


0.86 


1.72 


4 


0.91 ±0.06 


0.84 ±0.17 


2.28 ±0.14 


4.48 ± 0.50 


1.06 


0.85 


1.02 


1.25 


1.11 


3 


0.96 ± 0.09 


0.66 ± 0.25 


3.34 ± 0.22 


6.11 ±0.59 


0.92 


0.70 


0.87 


1.64 


1.02 


2.5 


1.21 ±0.14 - 


-0.12 ±0.39 


5.38 ± 0.42 


8.75 ± 0.73 


0.81 


0.44 


0.73 


2.18 


0.95 


2 


2.73 ±0.18 - 


-6.44 ±0.74 


21.0 ±2.3 


27.3 ±2.7 


1.34 


6.19 


2.36 


4.14 


2.72 
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Table 12: Fits for E$ with free parameters D u , Dd and C u , Cd subject to the constraint (|65|) . In all 
fits Cd takes it maximum allowed value 2.59 GeV -2 . Further explanations on the fits for E v are given 
in the text. 







fit parameters 






X 2 /data points 




X 2 /d.o.f. 


Pu 


Pd 


D u 


D d 


C u 


low 


high 


P 


n 


total 


4 


5 


0.38 ±0.12 


-0.65 ±0.13 


1.22 


1.89 


0.89 


1.28 


0.83 


1.32 


5 


5 


0.14 ±0.07 


-0.65 ±0.05 


0.85 ±0.14 


2.03 


0.90 


1.34 


0.79 


1.36 


6 


5 


-0.06 ± 0.06 


-0.65 ±0.05 


0.60 ±0.13 


2.13 


0.94 


1.40 


0.76 


1.41 


4 


6 


0.38 ± 0.07 


-0.82 ±0.04 


1.22 


1.94 


0.87 


1.29 


0.79 


1.31 


5 


6 


0.14 ±0.07 


-0.82 ±0.04 


0.82 ±0.14 


2.03 


0.91 


1.35 


0.76 


1.36 


6 


6 


-0.05 ±0.10 


-0.82 ±0.12 


0.58 ±0.46 


2.14 


0.94 


1.41 


0.74 


1.41 



Table 13: Fit for E v with fixed parameters C u = 1.22 GeV" 2 , C d = 2.59 GeV~ 2 and p d - p u = 1.60. 



fit parameters 




X 2 /data points 


X 2 /d.o.f. 


Pu D u 


D d 


low high p 


n total 


3.99 ±0.22 0.38 ±0.11 


-0.75 ± 0.05 


1.92 0.88 1.29 


0.80 1.31 



Table 14: Four-parameter fits for E$ with equal values of the parameters for u and d quarks. In the 
fit of the first row the constraint C q < 1.22 GeV -2 from positivity is imposed, whereas in the fit of 
the second row C q is left completely free. 



fit parameters 






X 2 /data points 




X 2 /d.o.f. 


a Pu = Pd 


D u = D d 


C u — c d 


low high p 


n 


total 


0.53 ±0.02 7.56 ±1.95 


-0.59 ±0.31 


1.22 


4.48 1.01 2.37 


13.26 


4.30 


0.61 ±0.05 2.24 ±0.71 


0.04 ±0.17 


3.50 ± 1.04 


3.87 0.80 2.00 


13.27 


3.97 



Table 15: Fits for E% with free parameters D u = Dd and C u , Cd subject to the constraint (1651) . 







fit parameters 






X 2 /data points 




X 2 /d.o.f. 


Pu 


Pd 


D u = D d 


C'U 


c d 


low 


high 


p 


n 


total 


7 


5 


-0.41 ±0.02 


1.22 


1.61 ±0.25 


2.79 


0.84 


1.61 


3.38 


1.98 


8 


5 


-0.49 ±0.02 


1.00 ±0.25 


1.78 ±0.28 


2.63 


0.91 


1.58 


1.88 


1.72 


9 


5 


-0.56 ± 0.02 


0.60 ± 0.20 


2.16 ±0.31 


2.58 


0.96 


1.60 


1.15 


1.62 


10 


5 


-0.63 ±0.02 


0.29 ±0.16 


2.48 ±0.32 


2.55 


1.06 


1.64 


0.76 


1.60 
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Table 16: Fits for with free parameters f3 u , (3d, D u , Dd and C u = Cd- In all fits C q takes its 
maximum value 1.22 GeV~ 2 allowed by the constraint 1)65(1 . 



fit parameters 




X 2 /data points 




X 2 /d.o.f. 




D u D d 


low 


high 


V 


n 


total 


8.97 ±2.86 6.75 ±1.82 


-0.62 ±0.34 -0.72 ±0.27 


2.34 


1.36 


1.75 


2.93 


2.12 


9.54 ± 1.79 6.49 ±1.27 


-0.69 ±0.20 


2.38 


1.35 


1.75 


2.93 


2.08 


7.76 ±1.23 


-0.47 ±0.18 -0.85 ±0.16 


2.24 


1.47 


1.77 


3.15 


2.13 


8.42 ± 1.30 


-0.74 ±0.16 


3.57 


1.11 


2.08 


15.9 


4.36 
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